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DMITRY CHELKAK A,B 

Abstract. We prove a number of uniform double-sided estimates relating discrete 
counterparts of several classical conformal invariants of a quadrilateral: cross-ratios, 
extremal lengths and random walk partition functions. These estimates hold true for 
any simply connected discrete domain (subset of a "roughly regular" planar graph, 
e.g., the standard square grid) with four marked boundary vertices, and are com- 
pletely independent of the domain geometry which can be very rough, having many 
fiords and bottlenecks of various widths. This allows one to use classical methods of 
geometric complex analysis for discrete domains "staying on the microscopic level" . 
Applications include a discrete version of the classical Ahlfors-Beurling-Carleman es- 
timate and some "surgery technique" developed for discrete quadrilaterals. 



1. Introduction 

1.1. Motivation. This paper was originally motivated by the recent activity devoted 
to the analysis of interfaces arising in the critical 2D lattice models (e.g., see [Smi06, 
SmilO] and references therein), particularly the random cluster representation of the 
Ising model [Kem09, DHN11, KS12, CDH12, CDHKS12]). At the same time, it has an 
independent interest, being devoted to one of the central objects of discrete potential 
theory on a (weighted) graph T embedded into a complex plane: partition functions of 
the underlying random walk running in a discrete (simply connected) domain Q. 

Dealing with some 2D lattice model and its scaling limit (an archetypical example 
is the Brownian motion in Q, which can be realized as the limit of simple random walks 
on refining square grids 5Z 2 ), one usually works in the context when the lattice mesh 5 
tends to zero. Then, it can be argued that a discrete lattice model is sufficiently close 
to the continuous one, if 5 is small enough: e.g., random walks hitting probabilities 
(discrete harmonic measures) converge to those of the Brownian motion (continuous 
harmonic measure, cf. [Kak44]) as 5 — > 0. After rescaling the underlying grid by 
statements of that sort provide an information about properties of the random walk 
running in large discrete domains Q C Z 2 . 

Unfortunately, this setup is not sufficient when we are interested in fine geometric 
properties of 2D lattice models (e.g., full collection of interfaces in the random cluster 
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representation of the critical Ising model): sometimes it turns out that one needs to con- 
sider not only macroscopic fi's but also their subdomains "on all scales (like 5 £ or even 
several lattice steps) simultaneously" in order to gain some macroscopic information. 
Questions of that kind are still tractable by classical means if those microscopic parts 
of Q are regular enough (e.g., rectangular-type subsets of Z 2 , cf. [KS12]). Nevertheless, 
if no such regularity assumptions can be made using some monotonicity features of the 
particular model, the situation immediately becomes much more annoying. 

Having in mind the classical geometric complex analysis as a guideline, in our 
paper we construct its discrete version "staying on a microscopic level" (i.e., without 
any passage to the limit or coupling arguments), which allows one to handle (possibly, 
extremely rough) discrete domains by more-or-less the same methods as continuous 
ones. We prove a number of uniform estimates (a "toolbox") which hold true for any 
simply connected Q, possibly having many fiords and bottlenecks of various widths, 
including very thin (several lattice steps) ones. Being interested in estimates rather 
than convergence, we do not need any nice "complex structure" on the underlying 
weighted planar graph. Instead, we assume that T and the corresponding random walk 
on it satisfy several rough "elliptic type" estimates (see Sect. 2.1 and 2.3 for details), 
deriving all results from those assumptions. 

It seems worthwhile to consider such "roughly regular" weighted graphs T, since 
we hope that our results can be applied in various setups (e.g., coming from lattice mod- 
els defined on nonregular lattices, circle packing constructions, or DLA-type processes 
in nonconstant environments) where some connection to discrete harmonic measure can 
be established (or is already plugged into the model). It should be said that we essen- 
tially use "uniformly bounded degree" and "quantitative local finiteness" assumptions 
(as well as uniform boundedness of edge weights), especially when proving a duality 
estimate for discrete extremal lengths. Thus, some important setups (notably, circle 
packings of random planar maps) are not covered, requiring some additional input 
(possibly, a kind of a "surgery" near high degree vertices, cf. [GN12]). On the other 
hand, results of our paper seems to be new even if T = Z 2 , so the reader not interested 
in full generality may always think about this, probably simplest possible, case. 

In order to shorten the presentation, below we widely use the following notation: 
assuming that all "structural parameters" of the graph T listed in Sect. 2.1 and 2.3 are 
fixed once forever (or if we work with some concrete planar graph, e.g., T = Z 2 ), 

• by " const " we denote positive constants (like ^ or 30 383 ) which do not depend 
on geometric properties (the shape of Q, positions of boundary points, etc) of 
the configuration under consideration or additional parameters we deal with 
(thus, " / ^ const " means that there exists a positive constant C such that 
the inequality / ^ C holds true uniformly over all possible configurations); 

• we write " / X g " if there exist two positive constants such that one has 
C±f ^ g ^ Cij uniformly over all possible configurations (in other words, / 
and g are comparable up to some absolute constants which we do not specify); 

• we write, e.g., "if / ^ const, then g 1 X </ 2 " iff, for any given constant c > 0, 
the estimate / ^ c implies C\g\ ^ gi ^ C^gi, where C\p. = C\^{c) > may 
depend on c but are independent of all other parameters involved. 
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1.2. Main results. The main objects of interest are (discrete) quadrilaterals, i.e., 
simply connected domains Q with four marked boundary points a, b, c, d listed coun- 
terclockwise. Focusing on quadrilaterals, we are motivated by two reasons. First, 
in the classical theory this is a "minimal" configuration which has a nontrivial con- 
formal invariant (e.g., all simply connected JTs with three marked boundary points 
are conformally equivalent due to the Riemann mapping theorem). Second, those are 
archetypical configuration for the 2D lattice models theory, where one often needs to 
estimate the probability of some crossing-type event between the opposite sides of Q. 

Note that, even if T = Z 2 , there is a crucial difference between discrete and contin- 
uous theories. The latter one is essentially based on conformal mappings and conformal 
invariance of various quantities, notably the conformal invariance of extremal lengths 
(see [Ahl73, Chapter 4] and [GM05, Chapter IV]). Using conformal invariance, one typ- 
ically may rewrite the question originally formulated in Q as the same question for some 
canonical domain (unit circle, half-plane, rectangle etc.), thus simplifying the problem 
drastically (e.g., see [GM05, Theorem IV. 5. 2]). In particular, up to conformal equiva- 
lence, (Q;a,b,c,d) can be described by a single real parameter (modulus). Therefore, 
all conformal invariants of those fTs (cross-ratios, extremal lengths, partition functions 
of the Brownian motion) are just some concrete functions of each other. 

This picture changes completely when coming down to the discrete level: for 
discrete domains (subsets of a fixed graph) we do not have any reasonable notion 
of conformal equivalence. Nevertheless, for a discrete quadrilateral, one can easily 
introduce natural analogues of all classical conformal invariants listed above. Namely, 
let = Zn([ab]n; [cd]n) denotes the total partition function of random walks running 
from the boundary arc [ab]n to the opposite arc [cd]n inside of f2 (if T = Z 2 , then Zq 
is the sum of weig hts 4- Lcngth W of all those nearest -neighbor paths 7, see Sect. 2.4). 
Then, we define the discrete cross-ratio Y n = Y^(a, b; c, d) of boundary points a, b, c, d 
as 

Y n := [Z n (a;d)Z n (b;c)/Z n (a;b)Z n (c;d)}^ 

(see Sect. 4.1 for details). We also use the classical definition of discrete extremal length 
(or, equivalently, effective resistance of the corresponding electrical network, see Sect. 6) 
Lfj = L^([a6]^; [cd]n) between [a6]n and [cd]a which goes back to Duffin [Duf62]. 

Certainly, one cannot hope that Zq,Yq and are related by the same identi- 
ties as in the classical theory. Nevertheless, one may wonder whether those can be 
replaced by some double-sided estimates which do not depend on geometric properties 
of (ft; a, b, c, d). One of the main results of our paper is Theorem 7.1 which gives a 
positive answer to this question. Namely, it says that, provided ^ const, one has 

ZflXYo and log(l + Y^ 1 ) x L n , 

uniformly over all possible discrete quadrilaterals. Note that we use discrete cross-ratio 
as an intermediary that allows us to relate "analytic" partition function and 
"geometric" extremal length in a way which is very similar to the classical setup. 

In order to illustrate a potential of the toolbox developed in our paper, we include 
two applications of a different kind. The first, given in Sect. 5, is a "surgery technique" 
for discrete quadrilaterals which is important for the fine analysis of interfaces in the 
critical Ising model, see [CDH12]. Namely, we show that it is always possible to cut 
Q along some family of slits L k into two parts Q' k and Q! k (containing [ab]n and [cd]n, 
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respectively) so that, for any k, one has 

Z n x Z n , k ([ab] n ; L k ) ■ Z n »(L fe ; [cd] n ) and Z n , k ({ab] n ; L k ) x k ■ Z n >>(L k ; [cd] n ) 

(see Theorem 5.1 for details). Using a discrete cross-ratios technique, we prove this, 
quite natural from a geometric point of view, result without any reference to the actual 
geometry of fl As always in our paper, double-sided estimates given above are uniform 
with respect to (Q; a, b, c, d) and k. 

Another application, given in Sect. 7, allows one to control a discrete harmonic 
measure Wdisc := ^n(u; [ab]n) of a "far" boundary arc [ab]n C dQ via an appropriate 
discrete extremal length L disc in f2 (see Sect. 7 and Theorem 7.8 for details). This should 
be considered as an analogue of the famous Ahlfors-Beurling-Carleman estimate (see 
[GM05, Theorem IV.5.2] and [GM05, p. 150] for historical notes). Again, we get a 
uniform double-sided bound which, as a byproduct, imply that 

log(l + Wdisc) X L disc ~ Lcont X log(l + W^) 

uniformly over all possible configurations (Q;u,a,b), where w C ont denotes the classical 
harmonic measure of the boundary arc [ab] seen from u in the polygonal representation 
of Q (see Corollary 7.9 for details). Note that results of this sort seem to be hardly 
available by any kind of coupling arguments. Indeed, dealing with thin fiords we are 
mostly focused on exponentially rare events for both discrete random walks and the 
(continuous) Brownian motion which are highly sensitive to widths of those fiords. 

1.3. Organization of the paper. In Section 2 we list several needed (rather mild) 
assumptions on the underlying weighted graph T and two important Assumptions (S) 
and (T) for the random walk on T (namely, uniform estimates for hitting probabilities 
and the expected exit time for discrete approximations of Euclidean discs). Then, 
we use (S) and (T) as a "black box", not discussing what one should ask about T to 
guarantee these properties. We also fix the notation for discrete domains Q, introduce 
the partition functions Zq of the simple random walk in Q, and discuss its relation 
with the standard notions of discrete harmonic measure and discrete Green's function. 

Section 3 is devoted to a uniform (up to absolute multiplicative constants) factor- 
ization of the three-point partition function Z^(a; [6c]q) via two-point functions Z^(a; b), 
Z^(a; c) and Zq(6; c). Namely, we prove that (see Theorem 3.5) 

Z n (a; [bc] Q ) x [Z Q (a; b)Z n (a; c) / Z Q (b; c) ]s 

uniformly over all configurations (Q; a, b, c). This is a cornerstone of our paper and the 
only one place where we involve some geometric considerations in the proofs. 

In Section 4, we introduce discrete cross-ratios Xq, Yq for a simply connected 
domain Q with four marked boundary points a, b, c, d (see Definition 4.3) and deduce 
from Theorem 3.5 several double-sided estimates relating X^, and Zq. In particular, 
we prove that X^ 1 x 1 + Y^ 1 (see Proposition 4.5, this is an analogue of the well known 
identity for classical cross-ratios) and Zn x log(l + Y^) (see Theorem 4.8, this is a 
precursor of the exponential-type estimate relating Zq and Lq). 

Section 5 is independent of the rest of the paper. It shows how one can use 
Theorem 3.5 and discrete cross-ratios introduced in Sect. 4 in order to build a sort of 
"surgery technique" which allows one to effectively "decouple" dependence Zq of the 
boundary arcs [ab]n and [cd]n by finding nice discrete cross-cuts in Q. 
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In Section 6, the notion of discrete extremal length Ln([a&]; [cd]) comes into 
play. We recall its definition and prove that is always uniformly comparable to its 
continuous counterpart - extremal length of the family of curves connecting [ab] and 
[cd] in the polygonal representation of Q (see Fig. IB). In particular, this fact implies 
the very important duality estimate for discrete extremal lengths (see Corollary 6.3). 
We also prove some simple inequalities relating Zq and L^ 1 (see Proposition 6.6). 

Section 7 summarizes all estimates for Yq, Zq and into the single Theo- 
rem 7.1 which is the culmination of our paper. Then, we show how to fit a discrete 
harmonic measure u>n(u; [ab]n) in this context (as Q \ {u} is not simply connected, a 
reduction similar to [GM05, p. 144] is needed). The main result (double-sided estimate 
of un(u; [ab]n) via an appropriate extremal length) is given by Theorem 7.8. As a 
simple byproduct, we prove Corollary 7.9 which says that the logarithm of a discrete 
harmonic measure is uniformly comparable to its continuous counterpart. 

Finally, in Appendix we derive several needed facts of discrete potential theory 
(elliptic Harnack principle, weak Beurling estimate, and pointwise estimates for Green's 
function) from Assumptions (S),(T). In some sense, we consider (S),(T) as a "pointe de 
la jonction" : being formulated for simplest possible domains (discrete approximations 
of Euclidean discs), they can be obtained by some specific means for various graphs, 
providing a basement for our toolbox which is more adapted for very rough JTs. 

Acknowledgments. The author would like to thank Hugo Duminil-Copin for several 
fruitful and motivating discussions, and many other colleagues, particularly Clement 
Hongler, Konstantin Izyurov, Antti Kemppainen, Eugenia Malinnikova, Asaf Nach- 
mias, Stanislav Smirnov and Mikhail Sodin for their helpful comments and remarks. 

This research was partly supported by P.Deligne's 2004 Balzan prize in Mathe- 
matics (associated research scholarship in 2009-2011), and the Chebyshev Laboratory 
at St. Petersburg State University under the Russian Federation Government grant 
11.G34.31.0026. Some parts of this paper were written at the IHES, Bures-sur-Yvette 
(October 2010), CRM, Bellaterra (May 2011), and Universite de Geneve (October 
2012). The author is grateful to these institutes for the hospitality. 

2. Notation, assumptions and preliminaries 

2.1. Graph notation and assumptions. Throughout this paper we work with an 
infinite undirected weighted planar graph (r; E r ) embedded into a complex plane C so 
that all its edges are straight segments (see Fig. 1) which is assumed to satisfy the list 
of assumptions (a)-(d) given below as well as Assumptions (S) and (T) formulated in 
Sect. 2.3. The notation T C C is fixed for the set of vertices which are understood as 
points in the complex plane (so \u — v\ means the Euclidean distance between u, v G T), 
and E r denote the corresponding set of edges. Each edge e G E r is equipped with the 
positive weight w e . Note that, in general, these weights are not related to the way how 
T is embedded into C. We assume that T satisfy 

• (a) uniformly bounded degree assumption: there is a constant vq > such that 
fj, v := J2(vv')eE r W 'W ^ V Q f° r all f G T and w e > v l for all e G E r . 

Clearly, this is equivalent to say that all edge weights w e are uniformly bounded away 
from and oo and all degrees of vertices of V are uniformly bounded as well. We then 
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denote random walk transition probabilities by 

W vv > w vv / 

T^vv' '■— = ^ • ^-J-J 

<L(W)gE r W vv> 

Note that zu vv > are uniformly bounded below by u^ 2 > and J2(w')eE r ZUm ' = 1, ^ G T. 
We now pass to the way how T is embedded into C. We assume that 

• (b) there are no "flat" angles: there exists a constant rjo > such that all angles 
between neighboring edges at any vertex v G T do not exceed ir — r]o (this also 
implies that all degrees of faces of T are uniformly bounded by 27r/?7o); 

• (c) edge lengths are locally comparable: there exists a constant po ^ 1 such 
that, for any vertex v G T, one has 

max \v ' — v I ^ Por v , where r v := min \v' — v\ (2.2) 

(™')eE r ' (^')eE r 

(below we sometimes call r v the local scale size); 

• (d) T is "quantitatively locally finite" : for any p ^ 1 there exists some constant 
u(p) > such that #{v G T : \v — u\ ^ pr u } ^ v{p), uniformly over all u G V . 

Roughly speaking, the last assumption prevents us from embedding Y so that vertices 
almost accumulate to a finite point, giving a uniform bound of their number in discs of 
radii comparable to the local scale size. Note that we do not assume that v(p) grows 
quadratically in p: the local scale size can vary from place to place but "not too fast" . 

Remark 2.1. It is easy to see that there exists a constant e > (depending on 
constants from (a)-(d) only) such that, for any u, v G T, the length ^22=1 \ u k+i~ u k\ of 
the shortest discrete path {uqU\. ..«„), (ukUk+i) G E r , connecting u = uq and v — u n is 
bounded by e$ 1 • \u—v\ (to construct a short path, use (b) to do the next step from uu 
towards v, if r Uk <^ \u k — v\, and (d), when reaching such that r Uk x \u k — v\). 

2.2. Bounded discrete domains. We start with a definition of a (bounded) discrete 
domain Q (see Fig. 1). Let (V n ; E^ t ) be a bounded connected subgraph of (r; E r ). In 
order to make the presentation simpler and do not overload the notation, we always 
assume that (ab) G E^ t for any two neighboring (in T) vertices a,b G V n (one can 
easily remove this assumption, if necessary). Denote by E^ d the set of all oriented 
edges (a int a) E^ t such that a int G V n (and a V n ). We set Q := Int Q U dQ, where 

IntQ := V n , dn := {(a; (a int a)) : (a int a) G E b n d }. 

Formally, the boundary dQ of a discrete domain f2 should be treated as the set of 
oriented edges (ai nt a), but we usually identify it with the set of corresponding vertices 
a, and think about Int fl and dQ as subsets of T, if no confusion arises. 

For a given vertex u G T and r > r u , we denote by B^(ti) the discrete disc 
of radius r around -u. Namely, IntB^(w) is the set of all vertices v G T lying in the 
connected component of T n {i> : \v — u\ < r} containing u (e.g., IntB^(w) = {«}), 
and dB>r(u) is the set of their neighbors, see Fig. 1A). 

Remark 2.2. It is worth to mention two useful corollaries of assumptions (a)-(d): 

(i) There exists a constant Eq > such that, for any u G V and r ^ r u , one has 
r fl {v : |i> — w| ^ £o r } C IntB^(-u) (this immediately follows from Remark 2.1). 
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Figure 1. (A) A discrete disc B^(u) and its polygonal representation 
(dashed). The inner vertices are colored black, the boundary ones are 
white. For all boundary edges (x^x) on the counterclockwise boundary 
arc [a6] B r( u ) C dB^(u) the corresponding midpoints x mi d are marked. 
(B) An example of a discrete simply connected domain Q with two 
marked boundary arcs [ab]n and [cd]^. Inside of Q, the notation used in 
the proof of Propostion 6.2 is shown: the polygonal representations of 
the neighborhoods A e and Aj of an edge and a face. 

(ii) For any given 5 > 0, there exists a positive constant v(25~ 1 ) > such that, for all 
u G T and r ^ r u , either r v ^ 5r for all vertices v G IntB^(w), or #B^(tt) ^ v{25~ v ) 
(this is a direct consequence of (d) as B^(-u) C B\ r {v)). 

We say that a discrete domain Vt is simply connected, if for any cycle in E^ t , 
all edges of T surrounded by this cycle also belong to Ef^ t . If £1 is simply connected, 
then its boundary vertices (or, more precisely, boundary edges) are naturally cyclicly 
ordered, exactly as in the continuous setting. For two boundary vertices a, b G dfl of a 
simply connected fl we denote a boundary arc [a6]n C dfl as the set of all boundary 
vertices lying between a and b (including those two) when one goes along dfl in the 
counterclockwise direction, see Fig. 1A (so, [a&^Uf&ajn = dfl and [a6]Qn[6a]^ = {a, b}). 
We also use the notation [ab)n := [ab]n \ {&}, (ab]n := [ab] n \ {a} etc. 

For a real function H : fl — >■ M. defined in fl, we define its discrete Laplacian by 

[AH](v) := E {vv ')eEr^w(H(v')-H(v)), v G Intfl, 

where the sum is taken over all neighbors of v, and w vv > are given by (2.1). We say 
that H is discrete harmonic in fl, if [AH](v) = at all inner vertices v G Intf2. Note 
that, if H > and [AH](v) = 0, then H(v')/H(v) ^ w^, ^ ^ ( see Sect 2 .1). 
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Below we often use two basic notions of discrete potential theory. The first is 
the discrete harmonic measure uJn{u\ E) of a boundary set E C dVt seen from 
an (inner) point u G fi. It can be defined as the unique function which is discrete 
harmonic in Q and coincides with on dQ. Note that uq(u; E) has a simple 

probabilistic meaning: it is equal to the probability of the event that the random walk 
having transition probabilities (2.1) started at u firstly hits dVl on E. The second is 
the (positive) Green's function Gn(v; u). It is the unique function which is discrete 
harmonic everywhere in Q except at u, vanishes on the boundary dQ, and such that 

[AG n (•;«)](«) = -|i" 1 . 

From the probabilistic point of view, Gq(v; u) is equal to the expected number of visits 
at u (divided by ji u ) of the random walk started at v and stopped when reaching dVL. 
Note that Gn is symmetric, i.e., Gn(u;v) = G^(v;u) (e.g., see Remark 2.6(h)). 

2.3. Random walk assumptions. We base our paper on two important assumptions 
for the random walk on T having transition probabilities (2.1). Namely, below we 
assume that two following properties hold true: 

Assumption (S) ("space"). There exist two positive constants f]o,CQ > such that, 
uniformly over all u G V , r ^ r u and 6 G [0, 27r], one has 

w B r( u )(u; {a G dBr(v) : arg(a-w) G [9, e + (n-r] )]}) ^ c . 

In other words, we assume that there are no exceptional directions: the random walk 
started at the center of any discrete disc B^(tt) can exit this disc through any given 
boundary arc of the angle n— rj with probability uniformly bounded away from 0. 

Assumption (T) ("time"). There exist two positive constants cq,Cq > such that, 
uniformly over all u G T and r ^ r u , the following is fulfilled: 

Despite (T) is formulated in terms of discrete harmonic functions only (which do 
not encounter any particular time parametrization of the random walk), it is natural to 
mention the following interpretation: let us consider some time parametrization such 
that the (expected) time spent by the walk at a vertex v before it jumps is of order 
rl (recall that local scales r v can be quite different for different v% if r 3> r u ). Then 
we ask the expected time spent in a discrete disc B r (u) by the random walk started at 
u before it hits dB r (u) to be of order r 2 , uniformly over all discrete discs. Note that, 
since X^eintB r (u) r l is comparable to r 2 due to (b), (c) and simple geometric reasons, 
(T) would follow from the following "analytic" claim: Green's function GB r ( u )( v ] u ) is of 
order 1 for those v G B r (u) which are away from u and dB r (u) and the total contribution 
of v 's lying close to u (where G-Q r ( u ) blows up) is bounded, see also Lemma A. 5. 

Remark 2.3. From now on, we think about all constants uq, rjo, po etc. used in assump- 
tions (a)-(d) and (S),(T) as fixed once forever. Thus, below we say, e.g., "with some 
uniform constants const i and const2" meaning that const i ; 2 may, in general, depend on 
u , r/o, po etc., but are independent of all other parameters involved (like domain shape, 
location of boundary points or particular graph structure). 
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Remark 2.4. In this paper, we use Assumptions (S) and (T) as a "black box", not 
intending to discuss what one should ask for the graph T and its embedding into C in 
order to guarantee them. Nevertheless, note that the following natural question arises: 

whether assumptions (a)-(d) are already sufficient for (S) and (T) or not? 

It is worth to mention that, if all lengths of edges e G E r are uniformly comparable 
to each other, then (S) and (T) for the simple random walk on V can be derived from 
the Parabolic Harnack Inequality technique (e.g., see discussion in [Koz07, Sect. 2.1]). 
An easier example when (S) and (T) can be obtained due to nice "local approximation 
properties" is given by some special random walks on isoradial graphs, see [CS11]. 

2.4. Partition function of the random walk in a discrete domain f2. The 

following notation generalizes both discrete harmonic measure and Green's function. 

Definition 2.5. Let Q C T be a bounded discrete domain and x,y E Q. We denote 
by Za(x; y) the partition function of the random walk joining x and y inside Vt. 
Namely, 

Z n (x;y) := £ w( 7 ), w( 7 ) := 1 } ^ = ^ vj UkUk+l , (2.3) 

ieS n (x;y) llk=0 k=0 

where Sq(x; y) = {7 = (u ~ u x ~ • • • ~ tt n ( 7 )) : u = x; u x , . . . , w„_i G Int Q; u n = y} 
is the set of all nearest-neighbor paths connecting x and y inside Q. Further, for 
A,B<zVL, we set 

and by RWo(A; B) we denote a random nearest-neighbor path 7 chosen from the set 
Sa(A; B) := U xe A, y eBSn(x; y) with probabilities proportional to w(7). 

Remark 2.6. It is easy to see that 

(i) if u G Int Q and b G dfl, then Zq(u; b) = fx^u^u; b); 

(ii) if both u, v G Int Q, then Zq(v; u) = Gq(v; u). 

Proof, (i) Focusing on the first step of 7 G Sn(u; b) in (2.3), one immediately concludes 
that the function 

"Zq(u; b), u G Int f2, 
fi^tlu — b], u G dfl, 

is discrete harmonic in Q and coincides with /i f J" 1 w^( • ; b) on the boundary dVt. Thus, 

Zq(u; b) = H{u) = u>n(u; b) for all u G Int Q. 

(ii) As above, it immediately follows from (2.3) that the function 



H(u) 



H{v) 



Zq(v; u), v G Int Q, 
0, v G dQ, 

is discrete harmonic everywhere in Q except at u, and 

H ( u ) = tK X + E(«u')eB r ^uu'Hiu'), 

where the first term /i" 1 corresponds to the trivial trajectory consisting of a single 
point u. Thus, [AH](u) = —fi^ 1 and Zq(v;u) = H(v) = Gn(v;u) for all v G Intfi. □ 
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3. Factorization theorem for the three-point function Z n (a; [6c] n) 

The main result of this Section is Theorem 3.5. It deals with a simply connected 
discrete domain Q and three marked boundary points a, 6, c G dfl (no assumptions 
about actual geometry of (Q; a, 6, c) are used) and provides a uniform up-to-constant 
factorization of the three-point function Z^(a; [6c]n) via Z^(a; b), Z^(a; c) and Z^(6; c). 
Actually, our proof is based on a factorization of the latter two-point functions via 
some inner point u G Int Q which is "not too close" to any of the boundary arcs [ab]a, 
[6c] n and [co\q. Thus, our strategy to prove Theorem 3.5 can be described as follows: 

• prove that the ratio Z^(a; u)Zq(u; 6)/Zn(a; b) is uniformly comparable with the 
probability of the event that RW$j(a; &) passes "not very far" from u (namely, 
at distance less than | dist(w; dfl)) - this is done in Proposition 3.1; 

• prove that this probability is bounded below, if u is "not too close" to any of 
boundary arcs [ab]n and [ba]n - this is done in Lemma 3.2 and Proposition 3.3; 

• find an inner vertex u which is "not too close" to any of [o6]q, [6c] q and [ca]n 
- this is done in Lemma 3.4 - and factorize all Z^'s using this u. 

Below we use the following notation. For a discrete domain Q and u G Int Q, let 

dn(w) := dist(«; dQ) = min |u-6|, B n («) := B^ (u)/4 (w). (3.1) 

Recall that (3.1) means IntB n (w) = {v G T : \v—u\ < \ dist(u; <9f2)} C Intfi (or, more 
accurate, a connected component of this set) and 8Bq(u) C Q is the set of all vertices 
neighboring to IntB^(-u). We also generalize notation (2.3) in the following way: for a 
given subdomain U C Q and a random walk path 7 = (u ~ U\ ~ • • • ~ «n( 7 )), we set 

■M7) : = ES^lKelntC/]- 
Then, for A, B C f2, we define 

ZnlT^KAjB) := £,^^(7)^(7). 
It is worth to note that 

=E[T p (RW n (A;B))] 

is the expected time spent in t/ by a (properly parameterized) random walk RW^(A; 5). 

Proposition 3.1. Let f2 6e a simply connected discrete domain, a, 6 G <9f2 7 and 
u G Int f2. T/ien, i/ie following double-sided estimate is fulfilled: 

Z J^ffl3 x p ( RWn (a; 6) n Int B n (ll) + 0] , 
Z^(a;6) 

some absolute (i.e., independent of Q,a,b,u) constants. 

Proof. Recall that both functions Zq( • ; a) = /i~ 1 w^( • ; a) and Zq( • ; 6) = /i fe " 1 wn( • ; b) 

are discrete harmonic and positive inside Q. Therefore, Harnack's Principle (see 
Lemma A. 2) gives 

Z n (u; a)Z n (u; b) x — (3.2) 



Et)6 



t>eIntB s2 (V) r v 

\2 



and EueintB n (u) r ^ x (dn( M )) due to simple geometric reasons. 
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Joining two random walk paths j av (from a to v) and 7„b (from v to b), and taking 
into account w^t^) = /i v ■ w(j av )w(-f vb ) x w(7 a „)w(7„ b ), it is easy to see that 

r^Z n (v;a)Z n (v;b) x Z n [T Bn(u) ](a;b) (3.3) 
t>elnt Bn(u) 

(indeed, each of the vertices Uk G RWo(a; b) contributing to TB n («) can be chosen as v 
in order to split RWn(a; b) into two halves 7 a „ and •y v b). 

Further, let w denote the first vertex Uk G Int Bq(u) of RW^(a; b), if such a vertex 
exists. Since in the right-hand side of (3.3) we do not count those paths which don't 
intersect Bn(u), by splitting RWn(ct;&) into two halves at w, (3.3) can be rewritten as 

Z n [T Bn(u) ]{a;b) x ^ Z n \ Bn (u)(a; w) Z n [T Bn(u) }(w; b), (3.4) 

u;6lntBn(u) 

where a (generally, doubly connected) discrete domain Q' := Q \ Bq(u) should be 
understood so that Int Q' = Int f2\Int B^(-u). It immediately follows from our definition 
of Z[T] and Harnack's Principle applied to the discrete harmonic function Zq( ■ ; b) that 

Zn[T Bn(u )](w;fc) - rlZ n (w;v)Z Q (v;b) x ^ r l z n(w; v) ■ Z Q (w; b) (3.5) 

t)6lntB n (tt) uelntBn(u) 

(indeed, for each w fc contributing to T Bn (u)(RW^(w; 6)) = Ylk=o r u k M u k e I n tBn(w)], 
split the random path RWq(w; b) into two halves "f wv , 7„& at the point v = Uk and use 
the up-to-constant multiplicativity ^{^wvlvb) x w(7 tut ,)w(7„b) once more). Moreover, 
using Assumption (T), it is easy to conclude that 

^ r^Z^(w;w) = Y r l G n(v] w) x (d^(w)) 2 for any u> G IntB^(-u). 

^elntBn(u) fGlntBn(M) 

(see Lemma A. 5 and Remark A. 6). Combining this with (3.2)-(3.5), one obtains 

Z n (u;a)Z n (u;b) 2 Z n [T Bn(u) ](a;b) 

— — ~ {(inW) — — 

Zn(a;6) Z n (a;b) 

_ E we intB n ( M ) Zn\Bn( M )(a; w) Z n (w; b) 
Z n (a;6) 

Finally, the numerator can be rewritten as 

Y Zn\B n (u){a;w)Z Q {w;b) x ^ w (7) 

u>€lnt Bn(u) 7GSn(a;&): 7nlnt B £2 (m)^0 

(as above, denote by w the first vertex G IntB^(-u) of 7, if such a vertex exists). 
Thus, (3.6) is comparable to the probability of the event RWn(a; b) D Bq(m) 7^ 0. □ 

Below we use the shorter notation 

P?[B n (u)] := P [RW^fa; 6) n Int B n (u) ^ 0] 

for the right-hand side of (3.1). Also, for wed and r ^ r u , we denote by B^(-u) the 
connected component of QnB^(u) containing u (and call it r -neighborhood in Q ofu). 



(3.6) 
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Lemma 3.2. Let Q be a simply connected discrete domain, a, b G dQ, and u G Intf2. 
Let L u denote the shortest discrete path running from u to dfl, d := Length(L M ) (recall 
that d ^ e^d^u), see Remark 2.1), and a,b G" B^u). Then, 

P [RW,(a; b) n L u ^ 0] ^ const • P^[B n ( M )]. 

Proof. Let w = m , «i, G Q be denned inductively by the following rule: Uk+i G L u 
is the first vertex on L u which doesn't belong to Int B^Uk) (thus, each -u fe+1 G <9Bn(ii fc ) 
and, finally, x := w n G <9f2). Further, let M m be the first of those w^'s such that the rest 
of the path L u is contained in the |dn(w)-neighborhood (in fl) of the ending point x. 

Since for any v G L u the distance from v to dQ along L u is comparable to da(v), 
one concludes that d^(uk) (or, equivalently, distances along L u ) decreases exponentially 
as k grows. In particular, this implies m ^ const. Let L' u and L" denote the parts of 
L u from u to u rn and from -u m to dVt, respectively. Harnack's principle implies 

H{u) x H( Ul ) x • • • x if (« m ) 

and so H x if(w) everywhere on L^ (with some absolute constants), for each of two 
functions H = Zq( ■ ; a) and H = Z^( • ; b). Therefore, 

P[RWnM)nL^0] < f;P^[B,M]xf: Z " (Mfc z a)Z ^ ;&) 

fe=o fc=o n ^ a ' ' 

< const. Za( "-f"'" ;6) xtf[B nM ]. (37) 

Due to Remark 2.2, without loss of generality, we may assume that r v $C 5od^(w) 
for all f G Bj n ( u )(a;), where 5o is chosen small enough independently of Q and u 
(otherwise, L" contains only a uniformly bounded number of vertices and the proof 
can be easily finished similarly to (3.7)). Then, it can be easily derived from (S) that 

H(v) ^ const -H(u) for all v G B" n(u)/2 (x) 

(see Lemma A.l and Remark A. 3). Moreover, since both functions H vanish on dQ 
near u n , the weak Beurling estimate (see Lemma A. 4) provides us a uniform bound 

H(u k ) ^ const .(dn(u fc )/dn(w)) A) ^(w), k ^ m, 

where /3q > is some absolute constant. In particular, the values H(uk) decay expo- 
nentially when k grows (as da(uh) do), and we easily finish the proof repeating the 
argument (3.7) given above: 

p[RW,M)nL^0] < g WWi 

— ' /. v ((i: b) 



k=m 



Z n (u;a)Z n (u;b) b 
^ COI1St 7 ( n u\ X P ^ l B ^( u )\ 

(the terms decay exponentially as k grows, so the final bound do not depend on n). □ 

Proposition 3.3. Let Q be a simply connected discrete domain, a, b G dQ, u G Int Q, 

and a > be such that both cun(u; [ab]n), ujq(u; [ba\a) ^ a. Then, the uniform estimate 

F a sj b [B n (u)} ^ const(a) (3.8) 



ROBUST DISCRETE COMPLEX ANALYSIS: A TOOLBOX 



13 



holds true, with some const(cr) > independent of Q,a,b,u. 

Proof. For simplicity, let us re-scale the underlying graph T so that dci(u) = 1. The 
weak Beurling estimate (Lemma A. 4) and our assumption on the harmonic measures 
of the arcs [ab]n and [ba]n seen from u imply that there exists big R = R(o~) (below we 
usually omit a, assuming that it is fixed throughout the proof) such that 

u is connected to both [ab]n and [ba]n inside B^(-u). 

Indeed, if u is not connected to, say, [ba]n, and R is big enough, then the harmonic 
measure u^(u; [ba]n) should be too small (recall that we set dn(u) = 1). 

Due to Remark 2.2(h), without loss of generality, we may assume that r v ^ 5 for 
all v G B2 R (u), where 5 is small enough (otherwise, the number of vertices in B^(-u) 
is uniformly bounded, and (3.8) is guaranteed by the fact that any random walk path 
connecting a and b in f2 should visit B^(-u) due to topological reasons). 

Moreover, it also follows from the weak Beurling estimate that there exist small 
r = r(a) > such that 

u remains connected to [ba]n in Vt 2 R t 2r ■— ^2r( u ) \ \J x e{ab] n ^2r( x ) 
(and vice versa with [ba] n and [ab]n interchanged). 

Indeed, if Uxe[afe] n ^2r( x ) separates u from [ba]n in B^ R ('u), then there exists a nearest- 
neighbor path 7 separating u and [foa]n in Q such that diam7 $C 5r and dist(w; 7) ^ |, 
which forces ujq(u; [ba]n) to be too small (see the second part of Lemma A. 4). 

Now let L denote some discrete path running from u to [ba]n in fl2R,2r- We define 
a sequence of vertices u = u ,Ui, ... G L inductively by the following rule: is the 
first vertex on L which does not belong to [J s<k ^ Bn(u s ). Let u rn be the first of those 
Uk satisfying dn(u m ) < s^r. In particular, the length of the shortest discrete path L Um 
connecting u m to dVt is less than r, see Remark 2.1. Hence, this path necessarily ends 
at the boundary arc (ba) n as u m U x e[o6] n ^2r( x )- 

Similarly to the proof of Lemma 3.2, we denote the parts of L from u to u m and 
from u m to dVt by V and L", respectively. Note that m is uniformly bounded. Indeed, 
by Remark 2.2(i) and our construction of u k , one has \uk—u s \ ^ So • ^dn(u s ) ^ for 
all ^ s < k ^ m. As all u k lie inside B^w), one has m ^ const. Applying Harnack's 
principle and Proposition 3.1 similarly to (3.7), one arrives at the uniform estimate 

m 

¥[RW n (a;b) n L' ± 0] ^ ^ ^ b [B Q (u k )\ < const • PS b [B n («)]. 

fc=0 

Finally, we apply Lemma 3.2 for the vertex u m (recall that d = Length(L Um ) ^ r 
and u m B^a) U B^ib) due to our choice of r). It gives 

P [RW^a; b) n L Um ^ 0] ^ const • F% b [B n (u m )] x F% b [B n (u)\. 

Summarizing, for a nearest-neighbor path L'U L Um running from u to (ba)n, we have 

¥[RW n (a;b) n (L'UL U J ^ 0] < const • P£ 6 [B n (u)]. 

Clearly, one can repeat the same arguments for the other boundary arc (ab)^. We 
finish the proof by saying that, due to topological reasons, RW^(a; b) should cross at 
least one of those two paths (joining u with (ba)n and (ab)n, respectively). □ 
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The last ingredient of the proof of Theorem 3.5 is the following simple 

Lemma 3.4. There exist a constant a > such that, for any simply connected discrete 
domain Q and three boundary points a,b,c £ dfl listed counterclockwise, one can find 
an inner point u £ Int Q so that all 

u n (u; [ab] a ), u n (u; [6c] n ), ^n(u; [ca] n ) > a . (3.9) 

Proof. Recall that "no flat angles" assumption (see Sect. 2.1) guarantees that all de- 
grees of faces of T are uniformly bounded. Let 

Intfi [a6 ] := {u £ Intfi : Lo n (u; [ab]n) ^ a}. 

If a is chosen small enough (independently of (f2; a, 6, c)), then VL\ab] is connected, con- 
taining all the vertices of faces touching [ab]n (in particular, d£l[ a b] D [a6]n). Moreover, 
due to the maximum principle, Q[ a b] is simply connected. Let 

L[ba] ■= dQ [ab ] \ [ab} n , 

b + £ L^] denote the next vertex on dQ[ a b] after 6, and a~ £ L[ ba ] denote the vertex 
just before a, if going along dQ[ a b] counterclockwise (in other words, b + and a~ are the 
endpoints of L^a] viewed as a boundary arc of f2[ a &])- For x £ L^a], let rc int £ Intf2[ a &] 
denote the corresponding interior vertex. Then, 

un(x int ; [6c] n ) + wn(^int; [ca] n ) ^ const • w n (x; [6c] n U [ca] n ) 

^ const ■ (1 — wj)(x; [a&]n)) ^ const -(1 — cr) 

for all x £ L^]. Moreover, for any two consecutive vertices x, x' £ L^j, the corre- 
sponding vertices Xi nt and x- nt belong to the same face of T, so 

wn(aj(nt; [6c] n) x wn(x int ; [6c] n ) and Wn(a4t5 t ca ]n) - Wn(^int; [ca] n )- 

On the other hand, wn(6^ t ; [6c] n) ^ const and wn(a^ t ; [ca]n) ^ const due to the same 
argument (e.g., 6^ t shares a face with 6i nt ). Therefore, observing Ltpa] step by step, one 
can find x £ L[ &a ] such that both u}n(x int ; [6c] n) and Co»n(x in t; [ca]n) are bounded below 
by some constant (and cun(xi nt ; [a6]n) ^ o" as x [nt £ Intfir a w). □ 

Theorem 3.5. Let Q be a simply connected discrete domain and a, b, c £ dQ be listed 
counterclockwise. Then, the following double-sided estimate is fulfilled: 

i 

(3.10) 

with some absolute (i.e., independent of Q,a,b,c) constants. 

Proof. Due to Lemma 3.4, one can find an inner point u £ Int Q such that (3.9) holds 
true, where <To is independent of Q, a, b, c. Since, for any x £ [6c]q, one has 

u n (u; [ax] n ) > Un(u; [ab} n ) ^ a and u n (u; [xa} n ) > u n (u; [ca] n ) > a Q , 

Propositions 3.1 and 3.3 imply 

Z n (a; [bc]a) = J2 xe [bc] n Z ^( a 5 x ) x £*e[bc]n Z ^( M 5 a )M u ', x) 

= Z n (u; a)Z n (u; [6c] n) x Z Q (w; a), 



Z n (a; [6c] n) 



Z n (a; b) Z n (a; c) 
Z n (6;c) 
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where we have used Zq(u; [bc]n) x ojq(u; [6c] n) x 1, and 



Z n (a; b) Z n (a; c) 
Zn(b; c) 



Z n (u; a)Zn(u; b) ■ Z n (u; a)Z n (u; c) 
Zq(m; b)Z n (u; c) 



as well. □ 

4. Discrete cross-ratios 

The main purpose of this Section is to obtain a uniform double-sided estimate (4.4) 
relating discrete analogues of two conformal invariants defined for a simply connected 
discrete domain Q with four marked boundary points a, b, c, d: discrete cross-ratio 
Yn(a, 6; c, d) (see Definition 4.3) and the total partition function Z^([a6]^; [cd]n) of 
random walks connecting two opposite boundary arcs. Note that the cross-ratio 
changes to its reciprocal when replacing boundary arcs [a6]n and [cd]n by "dual" ones 
([6c]n and [da]n), while the corresponding change of Zq is more sophisticated (see (4.4)). 

4.1. Monotonicity of Z (a;; a)/Zn(x; b) on the boundary and definition of 
discrete cross-ratios. Let two points a, b (or, more generally, two disjoint arcs A = 
[aia2]n and B = [6162] n) 011 the boundary of a simply connected discrete domain Q be 
fixed. Then, one can use the ratio Zn(x; o)/Zq(x; b) in order to "track" the position of 
x with respect to a, b. Being considered on dQ, this ratio has a monotonicity property 
(see Lemma 4.1 below), which allows one to use it as a "parametrization" of <9f2 between 
A and B. Namely, for x e <9fi, denote 

Lemma 4.1. Let Q be a simply connected discrete domain and A = [aia 2 ]n, B = [6162]^ 
denote two disjoint boundary arcs of Q. Then, Rq( • ; A, B) decreases along the bound- 
ary arc [a 2 6i]n and increases along the boundary arc [6 2 ai]n. 

Remark 4.2. In particular, if A = {a} and B = {b} are just single boundary points, 
then Rf}( • ; a, b) attains its maximal and minimal values on dQ at a and 6, respectively, 
being monotone on both boundary arcs [a6]n and [6a] q. 

Proof. Similarly to the proof of Remark 2.6(i), for any given t > 0, we define a discrete 
harmonic (in Q) function 



HAu) 



Z n (u; A) - tZ n (u; B), uelntQ, 

ji-\\ A (u) -tl B (u)), uedtl. 

Note that, for any x G <9f2, one has 

Z n (x; A) - tZ n (x; B) = H t (x) + w xx . nt H t (x- mt ). 

For a given boundary point x G (0261]^, let t x > be chosen so that H tx (xi Qt ) = 
(if x G (a 2 6i)n, this means R n (x;A,B) = t x as H tx (x) = 0, while Ra(bi, A, B) < t bl ). 

The function H tx is discrete harmonic in Q, vanishes on dQ \ (A U B), is strictly 
positive on A and strictly negative on B. Therefore, there exists a nearest-neighbor 
path running from x int to A such that H tx ^ along 7^. Due to the maximum 
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principle, this implies H tx (y int ) ^ for all intermediate boundary points y G [a 2 x)n. 
In other words, 

Zn(y, A) - t x Z n (y; B) = /x" 1 1 [y = a 2 ] + zu yy . nt H tx (y int ) > for all y G [a 2 x) n . 

Thus, Rn(y; A, B) > t x > Rn(x; A, B) for all y G [a 2 x)n which means that Rq(- ;A,B) 
decreases along [a 2 &i]n- The proof for the other boundary arc [&2ai]n is similar. □ 

Definition 4.3. Let Q be a simply- connected discrete domain and boundary points 
a, b,c,d G dfl be listed counterclockwise. We define their discrete cross-ratios by 



X n (a,b; c,d) :~- 



Z n (a; c) ■ Z n (b; d) 
Z n (a; b) ■ Z n (c; d) 



Y n (a,b; c,d) :~- 



Z n (a; d) ■ Z n (b; c) 
Z n (a;b)-Z n (c; d) 



Remark 4.4. Since a,b,c,d are listed counterclockwise, Lemma 4.1 implies 



X n (a,6; c, d) 



Rn(a; c, b) 
Rn(d] c, b) 



^ 1 and 



Xn(a, b; c, d) 
Yo(a, b; c, d) 



Rn{a; c, d) 
Rn(b; c, d) 



€ 1. 



We include the exponent | in Definition 4.3 by two (clearly related) reasons: 
first, it simplifies several double-sided estimates given below, and, second, it makes 
the notation closer to the standard continuous setup. Indeed, the continuous analogue 
of the partition function Z^(a;6) for the upper half-plane EI (up to a multiplicative 
constant) is given by {b — a)~ 2 , so the quantities and introduced above are 
"discrete versions" (in fl) of the usual cross-ratios 

(b— a)(d— c) . . (b— a)(d— c) 



x m (a,b; c,d) :~- 



and y M (a,b;c,d) :~- 



[c—a)(d—b) (d—a)(c—b) 

In the continuous setup, the following is fulfilled: (ii(o, b; c, d))^ 1 = l+(ym(a>, b; c, d)) -1 . 
One clearly cannot hope that the same identity remains valid on the discrete level for 
all fl's (even, say, if T is the standard square grid). Nevertheless, below we prove that 
the similar uniform double-sided estimate holds true for the discrete cross-ratios, with 
constants, in general, depending on parameters fixed in assumptions (a)-(d),(S),(T), 
but not on the configuration (fl; a, b, c, d) or the underlying graph Y structure. 

Proposition 4.5. Let fl be a simply connected discrete domain and a, b,c,d G dfl be 
listed counterclockwise. Then the following double-sided estimate holds true: 

(Xn(a, b; c, d))' 1 x 1 + (Y n (a, b; c, d))"\ (4.1) 

with some absolute (i.e., independent of fl,a,b,c,d) constants. 

Proof. We apply factorization (3.10) to both sides of the trivial double-sided estimate 

Z n (a; [bd\n) x Z n (a; [bc] n ) + Z n (a; [cd] n ), 

which is almost an identity besides the term Z^(a; c) counted once in the left-hand side 
and twice in the right-hand side. Dividing by [Z^(a; 6)Z^(a; c)Z^(a; d)} 1 ^ 2 , one obtains 
the following double-sided estimate: 



_Z n (a; c)Z n (b; d) 
which is equivalent to (4.1) 



Z n (a; d)Z n (b; c) 



+ 



Z n (a; b)Z n (c; d) 



(4.2) 
□ 
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Remark 4.6. It immediately follows from (4.1) that X n (a,b;c,d) x Y n (a,b; c, d), if 
Y n ^ const (which means that arcs [ab]a and [cd]n are "not too close" in Q). Moreover, 
the next Proposition shows that, in this case, Zn([a6]o; [cd]n) x Y^(a, b;c, d) as well, 
since is always squeezed (up to multiplicative constants) by and Y$> 

Proposition 4.7. Let Q be a simply connected discrete domain and a,b,c,d G dfl be 
listed in the counterclockwise order. Then, the following estimates are fulfilled: 

const -X n (a, b; c, d) ^ Z^([ab} n ; [cd] n ) < const -Y n (a, 6; c, d), (4.3) 

witt some absolute (i.e., independent of Q,a,b,c,d) constants. 

Proof. Due to Theorem 3.5, one has 

Zn(Hn; [cd]n) = Z n (x; [cd] n ) ~ / — r (Z n (x; c))*(Z n (a;; d))* . 

It follows from Lemma 4.1 that, for any x G [ab]n, 

(7 i whew i A\\h ( z Ux;c))^ (Z n (a;c))2 

(Z n (a:; c)) a (Z n (a:; d)) * = —j ■ Z n {x; d) > — i • Z n {x; d) . 

(Z n (a;;d))2 (Zn(a;d))a 

Therefore, summing and applying Theorem 3.5 once more, one obtains 

Z n (Hn; [crf]n) ^ const • — ^fe^ • Zn([a6] Q ; d) 

(Z n (c; d))2(Z n (a; d))2 

x (Z n(a;c ))j(Z n ( M ))l 

(Zn(c-,<0)5(Zn(o;6))i 

On the other hand, Cauchy's inequality (and Theorem 3.5 again) gives 

(Z n ( aft n; [cd\ n )) ^ const — — 

Zn{c;d) 

(Z n (a; c)Z n (6; c)Z Q (a; d)Z Q (6; d))3 



Z^(c; d)Z Q (a; b) 
= X n (a, b; c, d)Y n (a, 6; c, d) ^ (Y n (a, 6; c, cf)) 2 . □ 

Theorem 4.8. Let Q be a simply connected discrete domain and a, 6, c, d G <9f2 fre foted 
counterclockwise. Then, the following double-sided estimate holds true: 

Z n ([a&] n ; [cd] n ) x log(l+Y n (a, 6; c, d)) , (4.4) 

toii/j some absolute (i.e., independent of Q,a,b,c,d) constants. 

Proof. Denote Yq := Y^(a, 6; c, d), := Xn(a, b; c, d), and let a constant M be chosen 
big enough (independently of (f2; a, 6, c, d)). If Y^ ^ M, Propositions 4.5,4.7 imply 

Z n ([a&]; [cd]) ^ const -X n x (l+Yn^Yn ^ (1 + M)" 1 - log(l+Y n ), 
Z n ([a&]; [cd]) ^ const -Y n ^ const • M[log(l + M)]- 1 • log(l+Y n ) 
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(with constants independent of M). Thus, without loss of generality, we can assume 
that Y n ^ M (i.e., [ab]a and [cd]n are "very close" to each other in fi). Let 



Rn(x) 



r> ( A\ ^nfac) r LI 

Rn(x; c, a) = — - — — , xe[ab\ n . 



Z n (x; d) ' 

Due to Lemma 4.1, Rq increases on [ab]n. Moreover, it follows from Proposition 4.5 
(or directly from (4.2)) that 

Z n (b; c)Z n (a; d) 



'Rn(b)~ 


i 

2 


Z n (b; c)Z n (a; d) 






_Z n (b; d)Z n (a;c) 



1 + 



Z n (a; b)Z n (c; d) 



1+Y 



Yr 



As any two consecutive boundary vertices x,x' G [a L ]^ belong to the same face of T, 
one has Z n (x; c) x Zq(x'; c), Zn(x; d) x Z n (x'; d) and 

1 ^ - , < const . 
Rn(a;) 

Therefore, provided that ^ M is big enough, one can find a number n x logY^ 
and a sequence of boundary points a = a ,ai, . . . ,a n = b such that 

Rn(ofc+i) 
4 ^ < const 

Rn(ofe) 

for all fc = 0, n — 1. This can be easily rewritten as 

i 

Rn(ak) 



const < 



Xn(afc, Ofc+i; c, a) < - , 



Rn(afe + i) 

or, due to Proposition 4.5, as Yn(afc, ak+i', c, d) x 1. Hence, if the constant M was 
chosen big enough, the estimate (4.5) implies 

Zn(Ka fc+ i]n; Mo) x 1 
for all = 0, ...,n — 1. This easily gives 

n-1 

Z n ([a6] n ; [cd] n ) x Z n ([a fe a fe+ i]n; [cd] n ) x n x log Y n . (4.6) 

k=0 

Combining estimates (4.5) for Y^ ^ M and (4.6) for Y^ ^ M, one arrives at (4.4). □ 

5. Surgery technique 

The main purpose of this Section is to illustrate how tools developed above can 
be used to construct cross-cuts of a simply connected discrete domain Q having some 
nice "separation" properties, without any reference to the actual geometry of Q. The 
main result is Theorem 5.1 which claims the existence of those "separators". In Propo- 
sition 5.2, we also give some simple monotonicity properties of such cross-cuts. 

More precisely, let A = [aia 2 ]n and B = [b^jn be two disjoint boundary arcs 
of a simply connected Q. We are interested in the following question: whether it is 
possible to cut Q along some cross-cut L into two simply connected parts Qa, &b, one 
containing A and the other containing B, so that 

Z n (A; B) x Z Ua (A; L)Zq b (L; B) ? (5.1) 
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Int Q%[k] := I u G Int SI : Z")^ > k 



Moreover, we are interested not only in a single cross-cut L, but rather a family L k = 
L^[k] such that, in addition to factorization (5.1), one has 

Zn A (A;L k )/Z nB (L k ;B)^k. (5.2) 

Note that both Z Qa (A; L fc ), Z nB (L k ; B) > Z n (A; B). Thus, if Z n (A; B) > 1, then 
(5.1) certainly fails. Due to the same reason, one cannot hope for (5.2), if k <C Zn(A; B) 
or k ^> (Zq(A; B)) -1 . However, being motivated by the continuous setup, one certainly 
hopes for the positive answer in all other situations and, indeed, Theorem 5.1 claims 
the existence of a "separator" [k] and provides a natural construction of this slit for 
any given Q, A, B and k. 

Namely, let discrete domains fif [fc] and VL^(k~ l ) be defined by 

z n («;£) 

Int^AT 1 ) := Lelntfi: > k' 1 } 

(we use square and round brackets to abbreviate ^ and > inequalities, respectively). 
Below we always work with fc's which are not extremely big or extremely small, so 
that IntiTff/c] contains all vertices of faces touching A, while Int fi^A; -1 ) contains 
all vertices near B. Then, both Q^[k] and ^(k^ 1 ) are simply connected due to the 
maximum principle applied to the discrete harmonic function Zq( ■ ; A) — kZn( ■ ; B). 
Further, we denote the set of edges 

L B A [k] = L A B {k- 1 ) := {(u A u B ) G Eg t : u A G Int^[fc], u B G Int^AT 1 )} 

(see Fig. 2A). According to our conventions concerning the boundary of a discrete 
domain, this set can be interpreted as a part of dQ^[k], as well as a part of dfl^(k^ 1 ). 

Theorem 5.1. LetVt be a simply connected discrete domain, A,B C dQ be two disjoint 
boundary arcs, Z := Zq(A;B), and k > be chosen so that both Q A '■= ^a[^1 an d 
fl B := Vl^ik^ 1 ) are connected (i.e., VLa contains all inner vertices around A while fl B 
contains all inner vertices around B). Then, 

(i) for any fixed (big) constant K ^ 1, the following is fulfilled: 

ifZ^K and K^ 1 ^ k ^ K , then the cross-cut L k := L^[fc] satisfies both 
conditions (5.1), (5.2), with constants depending on K but independent 
of Q, A, B and k; 

(ii) there exists a (small) constant kq > such that the following is fulfilled: 

if Z ^ Kq and k^Z $C k ^ KqZ^ 1 , then the cross-cut L k satisfies both 
conditions (5.1), (5.2) with some absolute constants. Moreover, in this 
case, both Qa and Q B are always connected. 

Proof. Since Zq(ua] • ) x Zn(u B ; • ), it is clear that 
Z n (u A ] A) Z n (u B ; A) 

^Tl m x ^Tl m x k for a11 u = ( u ^b) e L fe . 5.3 

Z n (u A ;B) Z n (u B ;B) 

Let 8VIa l~l dQ, = [yA%A\n and dVt B C\dVl = [x B y B ]a (see Fig. 2A), and let 

Z A := Z n (A; [x B y B ]n), Z B := Z n {B; [y A XA}vi) , (5.4) 
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"1 A 




(B) 

Figure 2. (A) A simply connected discrete domain split into two 
parts [A;] and VL^ik" 1 ) according to the ratio of harmonic measures 
of two marked boundary arcs A = [aia 2 ] and B = [b^]. All edges 
(uaUb) that cross the slit L^[k] are marked, as well as four bound- 
ary edges (x^x A ), (x^x b ), (v^Vb), (vTva) e dQ neighboring to Lf[k]. 
(B) Notation used in Proposition 5.2 and schematic drawing of the mono- 
tonicity property Q^N c ^A UC N c ^aM for x e ( a 2&i)- 



where these partition functions are considered in the original domain fl. Then, 
Z QA {A;L k ) = Zn A (A;u B ) x ^ Z Qa (A; u B )Z n (u B ; dO) 



u£L k 



ueL k 



= ^2 Zn A (A;u B )Z n (u B ; [x B y B ]n) + ^ Zn A (A;u B )Z n (u B ] [vaXa]®) , 

u£L k wGLfe 

since Zq(u b ; dfl) x 1 for any u B G Intf2. Note that the first term can be rewritten as 

E«eL fc Z QA (A;u B )Z n (u B ; [x B y B ] n ) x Z n (A; [x B y B ] n ) = Z A . 

Indeed, each random walk path running from A to [x B y B ]fi inside Q should pass through 
L/j by topological reasons, so, denoting by u the first crossing, one obtains the result. 
Similarly, the second term is comparable to the total partition functions of those ran- 
dom walks which start from A, cross Lj, (possibly many times), and finish back at 
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[i/a^a]^- Denoting by v the last crossing of L k and using (5.3), one obtains 
^2 Zn A (A;u B )Z n (u B ; [y A x A ] n ) x ^2 Zn(A;v B )Z nA (v B ; [y A x A ]a) 

ueL fc v€L k 

x k ^2 ^(B; v B )Z nA (v B ; [y A x A ] n ) ~ kZ B , 
veh k 

since each random walk path running from B to [yA^n inside Q should cross L k . 
Thus, we arrive at the double-sided estimates 

Z nA (A; L fc ) x Z A + kZ B , 

and, similarly, Zn B (Lk; B) x k~ x Z A + Z B . Therefore, it is sufficient to prove that 

Z A /Z B x k and Z A Z B x Z . (5.5) 

It directly follows from (5.3) that 

Z n (x;A) Z Q (y;A) 

Z n (x;B) ~ ~ Z n (y; B) 1 ' } 

(here and below we omit subscripts of x and y, all the claims hold true for both 
x = x A ,x B and, similarly, y = y A ,y B , since the values of Zq(x a ; •) and Zq(x b ; •) are 
uniformly comparable). Denote 

Ya ■= Y n (ai, a 2 ; x, y), Y B := Yn(h, b 2 ; y, x), 
X A := Xn(ai, a 2 ; x, y), X B := X n (&i, b 2 ; y, x), 

where all discrete cross-ratios are considered in the original domain Vt. Using Theo- 
rem 3.5 and (5.6), it is easy to check that 

i 

„. /I W_ („.. /1M 2 

x fc. (5.7) 

The rest of the proof is divided into three steps: 

• first, we prove (5.5) assuming that both Z A ,Z B are bounded above by some 
absolute constant (roughly speaking, this means that x and y are "not too 
close" to both A, B), in some sense this is the most conceptual step which is 
based on discrete cross-ratios technique developed in Sect. 4; 

• second, we use discrete cross-ratios technique once again to show that, indeed, 
one has Z A ,Z B ^ const, if k x 1 (in particular, this implies (i)); 

• last, we analyze general case in (ii) by starting with k — 1 and then increasing 
it until Z A becomes x 1, which, as we show, cannot happen until fc max x Z _1 . 

Step 1. The proof of (5.5) under assumption Za,Zb ^ const. In this case 
Theorem 4.8 guarantees that Y A ,Y B ^ const as well, and Remark 4.6 says that 

Z A x [Y A X A ] 1 ' 2 and Z B x [YbXb] 1 ' 2 . 

Therefore, (5.7) immediately gives the first part of (5.5). Moreover, one has x 
and x Y^, which is equivalent to say that 

Zn(x;ai) _ Z n (x; a 2 ) ^ Z n (x;bi) _ Z Q (x; b 2 ) ^ 
Zfi(y;ai) ~ Z n (y;a 2 ) Z n (j/;&i) ~ Z^(y;6 2 ) ' 



[Y A X A ~ 


i 

2 


~Z Q {x; A)Z n (y;A) 


Y B X B _ 




_Z n (x;B)Z n (y;B) 
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In addition, Theorem 3.5 applied to (5.6) gives 

Z n (x; ai)Z^(rr; a 2 ) _ Z n (x;b 1 )Z n (x;b 2 ) 
Z n (y; a 1 )Z n (y; a 2 ) ~ Z n (y; &i)Z n (j/; 62) ' 

thus upgrading (5.8) to 

Z^Qe; oi) _ Z n (a:; q 2 ) _ Z^Qr; h) _ Z^Qr; & 2 ) ^ g ^ 

Zn(y;ai) ~ Z n (y;a 2 ) ~ Z n (?/;6i) ~ Z n (y;b 2 ) ' 
As Z ^ const, we also have Z x X^(ai, a 2 ; &i, 6 2 ). Rearranging factors, one obtains 

Z A Z g ^ [YaXaY^X^]^ ^ R 1 R Z n (a j ;x)Z n (x;b j )Z n (b j ;y)Z n (y;a j ) 
Z X n (a 1 ,a 2 ;b 1 ,b 2 ) 1 2 ' ^ (Zn(%; bj)Z n (x; y)) 2 

Finally, it directly follows from (5.9) that Yn(aj, x; 6j, y) x 1. Due to Proposition 4.5, 
this also implies Xfi(a,j,x;bj,y) x 1 and, similarly, Xn(x,bj;y,aj) x 1. Therefore, 

#j = [Xn(aj,x;&j,y)X n (x,&j;y,aj)]~ 3 x 1, 
i.e., Z^Zb x Z (which is the second part of (5.5)) and we are done. 
Step 2. Proof of Z A , Z B < const, if k x 1. In this case, Proposition 4.5 and (5.7) 
give 

YKI+Ya)' 1 x Y A X A x Y B X B x YKl+Ys)- 1 . 

Thus, if, say, Y^ ^ const, then Y B ^ const as well, and Z A ,Z B ^ const due to 
Theorem 4.8. Hence, without loss of generality, we may assume that both Y A ,Y B are 
bounded away from zero, which is equivalent to say that both X^,X B x 1, i.e., 

Z n (x; ai)Z n (y; a 2 ) x Z n (x; y)Z n (a X ] a 2 ) , 
Z n (x; b 2 )Z n (y; b ± ) x Z n (x; y)Z a (b 1 ; b 2 ) . 
Using Theorem 3.5 and (5.6), we obtain 

Zn(x; a 2 ) _ Z n (x; a 2 )Z n (x; Oi) _ (Z n (x; A)) 2 
Z n (y;a 2 ) Z n (ai, a 2 )Z n (x; y) Z n (x;y) 

_ (Z n (x; B)) 2 _ Z n (x; &i)Z n (x;& 2 ) _ Z n (x; 61) 
Z n (x;y) ~ Z n (6 1 ;fo 2 )Z Q (a;;y) ~ Z n (j/;&i) ' 
Due to Proposition 4.5 applied for a discrete quadrilateral (Q; a 2 , x; y, b x ), it gives 
Z n (x; a 2 )Z n (y; b x ) x Z n (x; y)Z n (a 2 ; b x ) x Z n (x; &i)Z n (y; a 2 ) . 

Similarly, one has Z n (:r; a 1 )Z n (y; b 2 ) x Z n (a;; y)Z n (ai, b 2 ) x Zn(x; b 2 )Z n (y; a x ). Then, 
using Xa,Xb x 1 and rearranging factors, one arrives at 

vv -vyvy ~ z n( fl i; 6 2)Zn(a 2 ;6i) _ v , , , v 

m'b — ia^-aibXb — 7T~? w 71 — r^r — in{ai,a 2 ,Oi,o 2 ) . 

Z n (ai; a 2 )Z^(6 2 ; 61) 

As Z is bounded above, Theorem 4.8 ensures that Y^(ai, a 2 ; 61, b 2 ) ^ const. Taking 
into account Y A ,Y b ^ const, we arrive at Y A ,Y b x 1, and so Z A , Z# x 1. 
Step 3. Proof of the general case in (ii). Let Z A (k) and Z B (k) be defined by (5.4) 
for a given k. Note that Z A (k), Z B (k) are piecewise-constant left-continuous functions 
of k which jump no more than by some constant factor Vq > 1 (see assumption (a) in 
Sect. 2.1), when domain £l A [k] (and, simultaneously, fi^(/c -1 )) changes. 
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We will fix k at the end of the proof, but in any case it will be less than 1. Since 
Z ^ 1, Step 2 ensures that Z^(1),Z B (1) ^ Co for some absolute constant Co (actually, 
Z^(l) and Z^(l) are much smaller, being of order Z 1 / 2 ). Now let us start to increase the 
parameter k. Since £1 A [k'] C Q A [k] f° r k' > k, the partition function Z A (k) increases, 
while Zb(/c) decreases. Let 

A; max := max{£; > 1 : Z A (k) ^ Co} • 

Due to Step 1, there exists a positive constant c ^ 1 such that the following is fulfilled: 

c k < Z A (k)/Z B (k) ^ c^k and c Z ^ Z A (k)Z B (k) ^ c X Z 

for any A; G [1,/Cmax]- Moreover, one has Zyi(/c max ) ^ z^Co) since the function Z A { • ) 
cannot jump too much at the point /c max . Therefore, we obtain the estimate 

, ^ Zyl^max) . 9 (Z^(A: max )) 2 -4a2 2 y-1 

K-max ^ <--0 ' ^ 71 7 ^ c ' 7} ^ U Q So c O ' n 

Thus, provided that k ^ min{l, t / o~ 4 Co c o}> 0*) holds true for all k G [lj/toZ^ 1 ] (and 
similar arguments can be applied for G [k ~ 1 Z; 1]). 

Finally, for all vertices near A, one has Z n ( ■ ; A) ^ const and Z ^ const -Zq( ■ ; B). 
Thus, choosing k small enough (independently of Q, A, B), one ensures that r2^[/t Z _1 ] 
is connected (and so Q A [k] is connected for all k ^ /toZ -1 as well). □ 

Dealing with more involved configurations (e.g. simply connected discrete do- 
mains with many marked boundary points, in addition to Theorem 5.1, it is useful to 
have some information concerning mutual "topological" properties of cross-cuts sepa- 
rating A and B, corresponding to different pairs A, B. In order to shorten the notation 
below, for x G dQ \ (A U B), we set 

n B M r>.B \T) ( a n\i f Zn(u;A) Z n (x;A)\ 

n A [x]:=n A [R n (x;A,B)] = l [ uen: ^-^ > ^-^ j. 

Roughly speaking, [#] is the set of those well which are "not further in Q" from 
A comparing to B than a reference point x. Note that, since the function R n ( ■ ; A, B) 
is monotone on the boundary arcs (a2&i)n and (62^1)0 (see Lemma 4.1), Q A [x] also 
behaves in a monotone way when x runs along dQ \ (A U B). 

Proposition 5.2. Let Q be a simply connected discrete domain, disjoint boundary arcs 
A = [aia 2 ]n, B = [fci^Jn andC = [ciC2]n & e listed counterclockwise, and BUC = [6iC2]n 
(i.e., b 2 and ci are consecutive points of dVL, see Fig. 2B). Then, 

fl A [x] C Q^lx] C Q A [x] for any x G (a 2 &i)n, 

n%\y] c QT c [y\ c ^SM M ^2/ y e (6 2 ai) n . 

Proof. Let x G (a 2 &i)n (the case y G (c 2 ai)n is similar) and -u G Intf^M- We need to 
check that u G Int Q A uC [x] . By definition, 

MGlnt^N Z n (u;A)-Z Q (x;C) ^ Z Q (x; A) ■ Z Q (u; C) , 
uelntQ^lx] & Z n (u;A)-Z n {x;BUC) ^ Z n (x; A) ■ Z n (u; BUC). 
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Since Z n ( • ; BUC) = Zq( • ; -B) + Z n ( • ; C), it is sufficient to prove that, for any b e B, 

Zn(x; b) = Z n (x; b- mt ) Z n (x;A) 
Z n (u;b) Z n (u;b int ) ^ Z n (u; A) ' 

For uffl, denote 

._ J Z ^( M > A ) ' Z ^( ;r ' u ) _ Z n (x; A) ■ Z n (u; v), v G Int ft, 

1 n~ x i[v— x], v g 9ft, 

Suppose that, on the contrary, H(b int ) < for some b E B. Since the function if is 
harmonic everywhere in ft except u (where it is subharmonic), and vanishes on 9ft 
everywhere except x (where it is strictly positive), there exists a nearest-neighbor path 
7f, M running from b mt to u such that H < along 7& u . On the other hand, H(c\ nt ) ^ 
for at least one c G C (otherwise, summation along the arc C gives a contradiction 
with the first part of (5.10)). Hence, there exists a nearest-neighbor path 7 CX running 
from Ci n t to x such that H ^ along 7 ca; . Since these two paths cannot cross each 
other and ft is simply connected, 7 c:E should separate u and A. Then, the maximum 
principle implies H(a int ) > for any a G A. Summing along the arc A, one arrives at 
the inequality 

Zn(u;A)-Zn(x;A)>Zn(x;A)-Z n (u;A), (5.11) 

which is a contradiction. 

Now let u G ft^ uC [a;]. Arguing as above, in order to deduce u G ft^[x] from the 
second part of (5.10), it is sufficient to prove that, for all c G C, 

Zn(x; c) _ Z n (x; c int ) Z Q (x; A) 
Z n (u;c) Z n (u;c iQt ) ^ Z n (u;A)' 

Suppose, on the contrary, that H(c- mt ) > for some c G C. Then there exists a path 
7 CX running from Ci nt to x such that H > along 7 ca; . Now there are two cases. If 
•y cx separates u and A, then the maximum principle implies H{a- m ^) > for all a <E A, 
which leads to the same contradiction (5.11). But if j cx does not separate u and A, 
then it separates u and B. Therefore, H(b int ) > for all b G B, which directly gives 
u G ft^[x] by summation along 5. □ 

6. Extremal lengths 

In this Section we recall the notion of a discrete extremal length Lo([a6]n; [cd\n) 
between two opposite boundary arcs of a discrete simply connected domain ft (which 
is nothing but the resistance of the corresponding electrical network), firstly discussed 
by Duffin in [Duf62]. Note that Lq can be defined in two equivalent ways: (a) via some 
extremal problem (see Definition 6.1), and (b) via solution to a Dirichlet-Neumann 
boundary value problem (see Proposition 6.4 and Remark 6.5). We also refer the 
reader to [BV12], where this equivalence and a connection to the Uniform Spanning 
Tree model are discussed in more details. The most important feature of (a) is that 
it allows one to estimate "in geometric terms". In particular, we show that 
is uniformly comparable to its continuous counterpart - extremal length of the corre- 
sponding polygonal quadrilateral (see Proposition 6.2 and Corollary 6.3 for details). 
At the same time, approach (b) allows us to relate Lq to the random walk partition 
function Zq discussed above (see Proposition 6.6). Note that this connection is of 
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crucial importance for the next Section, which starts with the complete set of uniform 
double-sided estimates relating Y n , Z n and (see Theorem 7.1). 

Let O be a discrete domain and E n = Efl t U E^ d be the set of edges of f2. For a 
given function ("discrete metric") g : E n — > [0; +oo) we define the "g-area" of f2 by 

where w e denote weights of edges of T (see Section 2.1). Further, for a given subset 
7 C E n (e.g., a nearest-neighbor path running in Q), we define its "g- length" by 

Finally, for a family £ of lattice paths in fl, we set L g (£) := inf 7£ £ L g (j). 
Definition 6.1. The discrete extremal length of the family £ is given by 

L[£]:= sup I^f, (6.1) 

where the supremum is taken over all g's such that < A g (Q) < +oo. In particular, if 
Q is simply connected, a, b,c,d G dQ are listed counterclockwise, and b ^ c, d ^ a, then 
we define Lo([a6]o; [cd]^) as the extremal length of the family (f2; [ab]n •<->■ [cd]^) 
of all lattice paths connecting the boundary arcs [ab]n and [cd] n inside Q. 

Definition 6.1 easily allows one to estimate the extremal length Ln([ab]n; [cd]n) 
from below, since for this purpose it is sufficient to take any "discrete metric" g in Q 
and estimate A g (Q) and L g (Q; [ab]n -H- [cd]n) for this particular g. Note that the most 
natural way to give an upper bound is to use (some form of) duality (6.6). 

For a (simply connected) discrete domain Q C T we denote its polygonal rep- 
resentation as the open (simply connected) set f2 c C C bounded by the polyline 
x mid x mid-" x mid x mid passing through all middle points a;^ id := |(x fc + xf nt ) of boundary 
cd ges (xf Qt x k ) G dQ in their natural order (counterclockwise with respect to Q), see 
Fig. 1. For a, b G dQ, a ^ b, we denote by [ab]^ C dfl c the part of this polyline from 
a m id to 6 m id, viewed as a boundary arc of Q c . In case a = b, we slightly modify this 
definition, setting, say, [aa]^ := [§(a~ id + a mid );a mid ] U [a mid ; |(a+ id + a mid )]], where 
a T denote the boundary points of Q just before and next to a. 

Let Lq := Lq([ci&]q; [cd]£j) denote the classical extremal distance between the op- 
posite arcs of a topological quadrilateral (Q c ; a mi d, 6 m id, c m i d , rfmid) in the complex plane 
(e.g., see [Ahl73, Chapter 4] or [GM05, Chapter IV)]). Note that our Definition 6.1 
replicates the classical one, which says 

r i 2 

Ln(Mn;Mn) = sup ± rr 2 ^, (6.2) 

g:QC^[ 0;+ oo) JJq9 dxdy 

where the supremum is taken over all g such that < JJ Q g 2 dxdy < +oo and the 
infimum is over all curves connecting [o6]q and [crf]^ inside f2 c (see [Ahl73, GM05]). 

Proposition 6.2. Let Q be a simply connected discrete domain and a, b,c,d G dQ, 
b 7^ c, d 7^ a, be listed in the counterclockwise order. Then, 

Ln([ab] n ; [cd] n ) x Lg([a&]H; [cd] € n ) (6.3) 

with some absolute (i.e., independent of Q,a,b,c,d) constants. 
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Proof. Let L disc := L n ([ab] n ; [cd] n ) and L cont := L^([a6]^; [cd\%). We prove two one- 
sided estimates separately, taking a solution to either discrete (6.1) or continuous (6.2) 
extremal problem, and constructing some related metric for the other one, thus ob- 
taining a lower bound for the other (continuous or discrete) extremal length. 

(i) L cont ^ const ■ Ldisc- Let g™ ax , e G E n , be the extremal metric in (6.1). For a face 
/ of T (considered as a convex polygon in C), let Af C T be defined by saying that Int Af 
consists of all vertices incident to /, and Aj be the polygonal representation of Af. 
Further, for an edge e G E r separating two faces / and /', let Int A e := Int Af U Int Ay 
and Ag be the polygonal representation of A e (see Fig. IB). We set 

g(z) := g^r; 1 !^), z G fi c , 

where r e denotes the length of e. Since each point in Q c belongs to a uniformly bounded 
number of edge neighborhoods A|p (recall that degrees of faces and vertices of T are 
uniformly bounded), one has 

fdxdy x (g?n 2 rl 2 Area(A^ n Q c ) x £ w e {g^)\ (6.4) 

as r^ 2 Area(A|p n f2 c ) x 1 x w e (see our assumptions on T listed in Sect. 2.1). 

Now let 7 be any continuous curve crossing Q c from [ab]^ to [cd]£j, F 1 be the set 
of all (closed) faces touched by 7, and E 1 C i? be the set of all edges of Q incident to 
those faces. It is clear that E 1 contains a discrete nearest-neighbor path running from 
[ab]a to [cd\ci- Thus, it is sufficient to estimate J^gds (from below) via ^2 eeE1 g™ ax . 
Note that, for any / G F 7 , 

7 should cross the annulus type polygon Ay \ / a£ /eas£ once. 

Let 7/ denote this crossing (there is one exceptional situation: if, say, b and c are two 
consecutive boundary points and / is a boundary face between them, then 7 may not 
cross the annulus Aj\f, so we denote by 7/ be the corresponding crossing of A^ itself). 
As degrees of vertices and faces of T are uniformly bounded, each piece of 7 belongs 
to a bounded number of 7/. Since Length(7/) ^ const -r e for any e ~ / (all those r e 
are comparable to each other due to our assumptions), we arrive at 

/ gds ^ const ■ ^ / gds ^ const • ^ Length(7/)o™ ax r~ 1 > const • ^ g\ 

Together with (6.4), this allows us to conclude that 

[inf 7 / qrfg] 2 [inf 7 L gmax (E7)]2 

Lcont ^ —77 — p , j > const — ^ const • L disc . 

JJ n g z dxdy A g m^{\l) 

(ii) L disc ^ const -L cont . Let g max : f2 c — >■ R + be the extremal metric in (6.2). It is 
well known that g max (,2) = \4>'(z)\ where conformally maps f2 c onto the rectangle: 

4> : Q c ^- {z : < Re z < 1, < Imx < L~„J, 

-1 -1 ( 6 - 5 ) 

a^L cont , &^0, C4l, d^l + tL cont . 



max 
ie 
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We set 



# max ds, e G E n . 



Note that, for each nearest-neighbor discrete path 7 in f2, we have Ee€ 7 9e = J 7 g max ds, 
thus it is sufficient to estimate E e gn e w e9e (from above) via JJ (g max ) 2 dxdy. 

Let z e denote the mid-point of an inner edge e. As is a univalent holomorphic 
function (in Ag nfi ), all values |0'(z)| for z G e are uniformly comparable to each other 
(and comparable to all other values |0'(z)| for z near z e ), e.g., see [Ahl73, Chapter 5, 
Theorem 5-3] or [GM05, Chapter 1, Theorem 4.5]. In particular, this implies 

g\ x r 2 e \<P\z e )\ 2 < const • I [ \<P'\ 2 dxdy. 

The same holds true for boundary edges: if Q c has an inner angle 9 X G (770; 2n] at the 
boundary point x mi d G dfl, then <fi behaves like (z — x mid ) 7T / 9x near x (or [z — x niid ) 7T / 2ex , 
if x is one of the corners a, b, c, d), hence \(f>'\ blows up not faster than \z— £ m id| _1//2 (or 
I -2- ^mid| -3 ^ 4 ) respectively) when 2 approaches x mi( i, which means g e x r e |0'(x int )|. 
As each point in fl c belongs to a uniformly bounded number of A^ , we obtain 



y w e g 2 ^ const • / / \<p'\ 2 dxdy 

Therefore, 

T ^ [mf 7 E eg7 ^] 2 ^ [inf 7 / 7 ^ S ] 2 
disc ^ "V w~^~ ^ COI1St " ff (V^lV^ ^ C ° nSt ' c ° nt • 



Corollary 6.3. Let Q be a simply connected discrete domain and a, b,c,d G dQ be four 
distinct boundary points listed counterclockwise. Then, 

Ln([ab] n ; [cd] n ) ■ L n ([bc] n ; [da] n ) x 1 (6.6) 

with some absolute (i.e., independent of Q,a,b,c,d) constants. 

Proof. Directly follows from (6.3) applied to both factors and the exact duality of 
continuous extremal lengths: L^([ab]^; [cd]£) • L 8(M8; [ da \n) = 1 - D 

We now pass to the second approach to the notion of extremal length via solution 
to the following Dirichlet-Neumann boundary value problem (which corresponds to the 
real part Re0 of the uniformization map (6.5)): 

Let Q be simply connected and a, b,c,d G dQ, b 7^ c, d 7^ a, be listed 
counterclockwise. Denote by V = V(n ; [a6]n,[cd]n) : ^ — > [0; 1] the unique 
discrete harmonic in Q function (electric potential) such that V = on 
[ab]n, V = 1 on [cd]n, and V satisfies Neumann boundary conditions 
(i.e., V{x iQt ) = V(x)) for x G dQ \ {[ab] n U [cd] n ). We also set 

(note that ^ xedn w XXiiit {V(x) -V(x int )) = E u6 intn A^[AV](u) = 0). 
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The next Proposition rephrase L^([a&]n; [cd]o) via I(V) (which is nothing but 
the electric current in the corresponding network). Note that, on the contrary to 
the classical setup, (6.7) does not allow to replace double-sided estimate (6.6) by the 
identity. Indeed, mimicking the continuous case, one can pass from V to its harmonic 
conjugate function V* that solves the similar boundary value problem for dual arcs, 
but this V* is defined on a dual graph T*, leading to the extremal length of some other 
discrete quadrilateral (drawn on T*) rather than ficT itself (see Remark 6.5). 

Proposition 6.4. For any simply connected discrete domain fl and any a, b,c,d G dfl, 
b 7^ c, d 7^ a, listed counterclockwise, the following is fulfilled: 

MHn; [cd] n ) = [I(V m ab] n ,icd\ (1 ))r [ - (6.7) 

Proof. See [Duf62] and [VorlO, BV12]. The core idea is to construct the function V 
using the extremal "discrete metric" g max for the family S := (Q; [ab]n -H- [cd]n) (which 
is unique up to a multiplicative constant due to simple convexity reasons). Namely, 
let (Q;u -H- [ab]ci) denote the family of all discrete paths running from u G VL to the 
boundary arc [ab]n inside Q, and 

V(u) := L ff max(Q;u -h- [a6]n). 

Then, V is constant on [cd]n and satisfies Neumann boundary conditions on (be) n, 
(da)ci (if one of these properties fails, then one can improve g max on the corresponding 
boundary edge so that L g {£) does not change while A g (Q) decreases). In particular, 
one can normalize g max so that V = 1 on [cd]n. 

Moreover, V is discrete harmonic in Q. Indeed, note that V(u') — V(u) = ±5^/ x 
for any {uu') G E Q (otherwise, one can improve g™*)- Then, for a given u G Int Q, 
replacing g™* by g™*+ £ 011 a ^ edges (uu') G E n such that V(u') > V(u) and, simulta- 
neously, replacing g™* by g™ x — e on all (uu') G E n such that V(u') <V(u), one does 
not change any global distances (and, in particular, does not change L g (S)), while the 
area A g (Q) changes by efj, u [AV](u) + 0(e 2 ). 

Finally, using discrete integration by parts and [AV](u) = 0, one concludes that 

L^ 1 = V(^) = J2e=(uu')eE^e(V(u')-V(u)) 2 

= -Euemtn^[AV}(u)V(u) - E x e d ^ i J V (^t)-V(x))V(x) 
= ^ [ ^x^(l-V{x int ))=I(V). □ 

Note that, for any discrete harmonic in Q function V, one can construct a discrete 
harmonic conjugate function V* which is uniquely defined (up to an additive constant) 
on faces of Q (including boundary ones) by saying 

H(0)-H(fT) ■= • (H(v')-H(v)) (6.8) 

for any oriented edge (vv') G E n , where and /^f ht denote faces to the left and to 
the right of (vv'), respectively. The function V* is well defined locally (iff Ay = 0), 
and hence well defined globally, as f2 is simply connected. Moreover, for any inner face 
/ in Q, it satisfies a discrete harmonicity condition 

E r ^fAv*(f)-v*(f)) = o, (6.9) 

where wjy := w~J, for any couple of dual edges (//') = (vv')*. 
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Remark 6.5. If one takes V = V(fi;[a6]n,[cd|n)j then the harmonic conjugate function V* 
is constant along boundary arcs (bc)n and (da)n (since V satisfies Neumann boundary 
conditions on these arcs). Fixing an additive constant so that V* = on (be) q and 
tracking the increment of V* along [ab]n, one obtains V* = I{V) on (da)n. Further, 
Dirichlet boundary conditions for V on [ab]n and [cd]n can be directly translated into 
Neumann ones for V* (one can easily see that V* satisfies (6.9) with smaller number 
of terms at all faces touching [ab]n or [cd]n). Thus, [/(V)] -1 • V* solves the same 
Dirichlet-Neumann boundary value problem for the dual quadrilateral drawn on T*. 
Moreover, 

E (//re ^ w //'(^(/0-^(/)) 2 = E {vv r } ^^AV(v')-V(v)r = La 1 , 

and hence the corresponding dual extremal length is equal to [/(V) -2 ]!^ 1 ] -1 = Lq 1 . 

The last Proposition in this Section gives an estimate for the partition function 
Zn((a&]n; [cd]n) of random walks joining [ab]n and [cd]n in Q via the extremal length 
Ln([a6]n; [cd]n) (note that the latter can be thought about as the (reciprocal of) similar 
partition function for random walks reflecting from the dual boundary arcs). 

Proposition 6.6. Let Q be a simply connected discrete domain, and a, b,c,d G dQ, 
b c, d a, be listed counterclockwise. Then, the following is fulfilled: 

Z n ([ab] n ; [cd]n) < const -(L n ([a6]n; [cd}n))~\ 

where const does not depend on Q,a,b,c,d. Moreover, ifLn([ab]n] [cd] n ) ^ const, then 

Z n ([a6] n ; [cd] n ) x (L n ([a6] n ; [cdja))' 1 

(with constants depending on the upper bound for Lq but independent of Q,a,b,c,d). 

Proof. It is easy to see that, for any u e Intfi, V(u) is equal to the probability of 
the event that the random walk started at u and reflecting from complementary arcs 
[bc)n,(da)n exists Q through [cd]n (indeed, this probability is a discrete harmonic 
function which satisfies the same boundary conditions as V). Hence, for any x e [ab]n, 

V(x int ) ^ const • Z n (x- mt ; [cd} n ) x Z n (x; [cd] n ), 

since the right-hand side is (up to constant) the same probability for the random walk 
with absorbing boundary conditions on (6c) n and (da)n- Thus, (6.7) gives 

(L n ([a6] n ; [cdjn)) -1 = J2 x e[ab] n w xx int V(x int ) > const • Z n ([ab] n ; [cd] n ). 

Further, let := L^Qafrj^; [cd]^) ^ const. Due to duality given by Corollary 6.3, 
it is equivalent to Lq([6c]q; [da\n) ^ const. We have seen above that this implies 

Z Q ([bc] Q ; [da]a) ^ const 

which is equivalent to Yn(b, c; d, a) < const due to Theorem 4.8. Therefore, 

Z n :=Z n ([ab];[cd]) x log(l+Y n (a, 6; c, rf)) = \og(l + (Y n (b, c; d, a))' 1 ) > const . 

Since Zq ^ const -L^ 1 in any case, this implies Zq x 1, if x 1. 

Thus, we are mostly interested in the situation when is very small (i.e., bound- 
ary arcs [ab]n and [cd]n are "very close" to each other in Q). Our strategy in this case is 
similar to the proof of Theorem 4.8: we split [ab]n into several smaller pieces [0^0^+1}^ 
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such that LnQafcafc+i]^; [cd\n) x 1 and apply the result obtained above to each of these 
smaller arcs. Recall that 

x£[ab] n 

We construct the boundary points a = ao, Qi, o n +i = b £ dQ inductively by the 
following procedure: if is already chosen, we move au+i further along the boundary 
arc [ab]n step by step until the first vertex a^+i such that 

(LQ([a fc a fe+ i]n; [crf]n)) _1 = w a . a!lIrt V r (n;[o fc o fc+ i]n 1 [cdin)(aJiiit) > 1 

(or afe + i = b). Note that this sum cannot increase by more than some absolute con- 
stant Vq on each step (as we increase the absorbing boundary [ofcdfc+ijn, all terms 
decreases, while the new (last) term is no greater than w xx . nt ^ u ). Therefore, 
LQ([afcdfc+i]fi; [cd]n) x 1 for all k, possibly except the last one (when we are forced 
to choose a n+ i = b before the sum becomes large). As we have seen above, this implies 

Zf}([a fe a fe+ i]n; [cd] n ) x 1 for all k = 0, 1, ...,n-l. 

Note that V = V { a ;[ab]n:[cd]n) ^ Vp. [a k a k+ i]n,[cd]n) due to monotonicity of boundary con- 
ditions (the absorbing boundary is larger in the first case), which gives 

^ W XXint V(x int ) ^ W xx int V(n;[a k a k+1 h,[cdh){Xint) < COnst 

for all fc = 0, 1, ...,n, thus L _1 ^ const -(n+1) x n, which implies the inverse estimate 

n 

Z n x ^ Zn([a fc a fc +i]n; [cd]n) ^ const -n x L^ 1 . □ 

k=0 

7. Double-sided estimates of harmonic measure. 

We start this Section by Theorem 7.1 which combines uniform estimates obtained 
above for cross-ratios Y^, partition functions Zq and extremal lengths of discrete 
quadrilaterals (Q; a, b, c, d). Then, we show how tools developed in our paper can be 
used to obtain exponential double-sided estimates in terms of appropriate extremal 
lengths for the discrete harmonic measure Uq(u; [ab]n) of a "far" boundary arc (similar 
to the classical ones due to Ahlfors, Beurling and going back to Carleman, see [Ahl73, 
§ 4-5,4-14] and [GM05, § IV.5,IV.6]). The main result is given by Theorem 7.8. In 
particular, it allows us to obtain a uniform double-sided estimate of \ogun(u; [ab]n) via 
\ogu)Qc(u; [ab]o), where Uqc denotes the continuous harmonic measure in a polygonal 
representation of Q (see Corollary 7.9). Note that one cannot hope to prove the similar 
estimate for 0Jn(u; [ab]n) itself: dealing with thin fiords, one faces with exponentially 
small harmonic measures which are highly sensitive to widths of those fiords. 

Theorem 7.1. Let Q be a simply connected discrete domain and distinct boundary 
points a, b,c,d G dQ be listed counterclockwise. Denote 

Y := Yn(a, b; c, d), Z := Z n ([ab] n ; [cd] n ) , L := L n ([ab] n ; [cd]n) , 
Y':=Y n (&,c;d,a), Z' := Z n ([6c] n ; [da] n ), L' := L n ([6c] n ; [da] n ). 
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(i) If at least one of the estimates 

Y ^ const, Z ^ const, L ^ const, 

Y' > const, 21 ^ const, L' < const ^ ' ' 

holds true, then all these estimates hold true (with constants depending on the initial 
bound but independent offl,a,b,c,d). Moreover, if at least one ofY, Y' , 2, 21 , L, L' is 
of order 1 (i.e., admits the double-sided estimate x lj, then all of them are of order 1. 

(ii) If (7.1) holds true, then the following double-sided estimates are fulfilled: 

ZxY and log(l+Y _1 ) x L. 
In particular, there exist some constants flip, > such that the uniform estimate 

d ■ exp[-/3iL] s$ Z sC C 2 • exp[-/3 2 L] (7.2) 
holds true for any discrete quadrilateral (fl;a,b,c,d) satisfying (7.1). 
Proof, (i) It follows from Theorem 4.8 and Proposition 6.6 that 

log(l+Y) x Z ^ const • L- 1 and log(l+Y') x Z' < const • (L') _1 . 
Moreover, YY' = 1 by definition, and LL' x 1 due to Corollary 6.3. Therefore, one has 

Y ^ const Z ^ const <= L ^ const 

t t 
Y' > const <^> Z' > const -<= L' < const, 

which gives the equivalence of all six bounds. Interchanging Y, Z, L and Y', Z', L', 
one obtains the same equivalence of inverse estimates. Thus, if at least one of these 
quantities is x 1, then all others are x 1 as well. 

(ii) Since Y ^ const, Remark 4.6 guarantees that ZxY. Further, as L' ^ const, 
Proposition 6.6 gives Z' x (L') _1 , and hence 

Ml+Y- 1 ) = log(l+Y') x 21 x (L')" 1 x L . 

Thus, exp[/3 2 L] ^ 1+Y" 1 ^ expfftL] for some (5 1>2 > 0, and 1+Y" 1 x Y" 1 x Z" 1 . □ 

Now let u G Int Vl and [ab]n C fi be some boundary arc of ft which should be 
thought about as lying "very far" from u (so that the harmonic measure 0Jn(u; [ab]n) 
is small). In order to be able to apply exponential estimate (7.2) to this harmonic 
measure, one should firstly compare the partition function of random walks running 
from u to [ab]n in ft with a partition function of random walks running between opposite 
sides of some quadrilateral. Recall that we denote by dn(u) the (Euclidean) distance 
from u to <9ft and = Bq(u) is a discrete disc around u of radius |dn(u): 

IntBn(u) := {v G Int ft : \v-u\ < \dn(u)}. 

Let a discrete domain = Aq(u) be defined by 

Int A n (u) := Int ft \ Int B n («). 

Note that, according to our conventions, one can identify the inner part of 8Aq which 
we denote by = Cq(u) (i.e., dAn = Cq U dft) with OBq. 
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(A) (B) 

Figure 3. (A) In order to analyze the discrete extremal length between 
Cq and [ab]ci, we cut a doubly connected domain along some nearest- 
neighbor path 7 running from c G to d E (ba)n, so that two identical 
copies of 7 are included into a simply connected domain A^ (which is 
drawn on the universal cover of An). Thus, the boundary <9A^ is formed 
by the outer part (cf ight d left ) = dQ, the inner part ( c lcft c right ) = 
and two paths 7 lcft and 7 n s ht consisting of vertices neighboring to 7. 
(B) If a vertex u is close to dfl, it might happen that A^ is simply 
connected or even disconnected. Then, we denote by A^ the proper 
connected component of A^, and by C' n the corresponding part of dA' n . 

Remark 7.2. Here and below (Lemma 7.3 - Theorem 7.8) we assume that Int Aq(u) is 
doubly connected (in order words, Bq(u) does not touch the boundary faces of dQ, i.e., 
u is not too close to dfl in the graph metric). Otherwise, one can apply an appropriate 
version of Lemma 7.3, which relates ojq(u; [ab]n) to the partition function of random 
walks running in Aq(m), and directly estimate the latter partition function by the 
corresponding discrete extremal length using (7.2) (e.g., see the proof of Corollary 7.9). 

Lemma 7.3. Let a simply connected discrete domain fl and u e Int Q be such that 
Aq(u) is doubly connected, and [ab]n C dQ. Then, 

u n (u; [ab] n ) x Z A ^ u) (C n (u); [ab] n ) ■ (7.3) 

Proof. For a random walk running from u to [ab]^ in Q, let v denote its last vertex 
on Cq (such a vertex exists due to topological reasons). Splitting this path into two 
halves (before v and after v, respectively), one concludes that 

u n (u; [ab] n ) x Z n (u; [ab] n ) x ^ Z n (u; v)Z An (v; [ab]n) 

veCn 

As Zq(u;v) = Gn(u;v) x 1 for any v G (see Lemma A. 5(h)), this gives (7.3). □ 

In order to relate the partition function (7.3) of random walks in the annulus 
Aq(u) to a partition function of random walks in some simply connected domain, 
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below we cut Aq(u) along appropriate nearest-neighbor paths 7 = (c int ~ ... ~ d int ) 
such that c <E Cn and d G dfl \ [ab]^. For a given 7 (which is always assumed to 
be a non-self-intersecting path on the universal cover of A^), we define a simply 
connected domain (see Fig. 3A) as follows: 

z/7 left ; 7 nght are two copies of 7 lying on consecutive sheets of A^, then 
Int A£ := 7 lcft U [(Int A n ) \ 7] U 7 risht C IntAg. 

In other words, we cut A^ along 7, accounting both sides of the slit as interior parts of 
a discrete domain A^ (which is, in particular, always connected and simply connected). 
We then denote by 7^ and 7^| the corresponding parts of <9A^, thus 

dAl = (d le V^) A o U 7^ U (c" gl V cft ) A o U 7 £f , 

where disjoint parts of <9A^ are listed counterclockwise with respect to A^. We also 
use simpler notation (<i left (i right ) A o = dVl and (c nght c lcft ) A o = Cq : if no confusion arises. 

Corollary 7.4. Let a simply connected discrete domain Q and u G Int Q be such that 
Afi(w) is doubly connected, and [ab]n C dQ. Then, for any nearest-neighbor path 7 
running from Cq(u) to (ba)n, the following is fulfilled: 

const • Z A 7 (C n ; [ab] n ) ^ u> n (u; [ab] n ) ^ const • Z A 7 (7^ UC.U 7 g ht ; [ab] n ). (7.4) 

Proof. Indeed, 

Zajt (C n ; [ab] a ) < Z As2 (C„; [o6] n ) ^ Z A , (7^ UC.U 7*?*; Ho) 

due to simple monotonicity properties of the random walk partition function Zq with 
respect to domain Q: e.g., for the left bound, one forbids the random walks running 
from Cq to [ab]n to cross 7 (still allowing them to touch 7 or to run along it). □ 

Theorem 7.1 (namely, (7.2)) allows one to estimate both sides of (7.4) using 
corresponding discrete extremal lengths. We now prove that one can choose 7 so that 
both those extremal lengths are comparable to the extremal length of nearest-neighbor 
paths connecting and [ab]n in the annulus An- 

Remark 7.5. Below we apply Propositions 6.2 and 6.4 to a doubly connected discrete 
domain A^ and its inner boundary instead of a second boundary arc [cd\n of a 
simply connected domain Q. It is worth to note that we didn't use any "topological" 
arguments in the proofs of those Propositions. 

Proposition 7.6. Let a simply connected discrete domain Q and u G Int Q be such 
that Aq(u) is doubly connected, and [ab]n C dQ. Then, 

(i) there exists a nearest-neighbor path 7 running from Cq to (ba)n such that 

L A -r(C n ;[a6] n ) ^ 2L An (C n ;[ab] n ); 

(ii) for any given q > 1, either L An (Cn; [a6]n) < Q 2 L As2 (C^; dfl) (i.e., the arc [ab]n is 
not so far from u), or there exists a nearest-neighbor path 7 running from to (ba)ci 
such that 

L A72 ( 7 J ) c d ft UCnU7bd ght ;[^]n) > (l-g-^LAjCnjMn). 
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Remark 7.7. (i) The constant 2 in the first estimate is an overkill: as it can be seen 
from the proof, both sides are almost equal to each other for a proper slit 7. 
(ii) Since discrete and continuous extremal lengths are uniformly comparable to each 
other, for any fl and u, one has 

L Af2 (Cn;^)xL Ag (C^;^ c )xl. 

Proof. Let V = V(A n ;[ab] n ,c n ) '■ — > [0; 1] be the unique discrete harmonic function 
such that V = on [a6]n, V = 1 on Cn, and V satisfies Neumann boundary conditions 
on dfl \ [ab]n. Recall that Proposition 6.4 (see also Remark 7.5) says 

(L An (Cn;[ab] n ))- l = I(V)= £ w XXint V (x int ) = £ w vv , (V (v') -V (v)) 2 . 

xe[ab]si (w')inA^ 

(i) Let V* denote a harmonic conjugate function to V (see (6. 8), (6. 9) and Remark 6.5) 
which is defined on the universal cover Aq of A^. Tracking its increment along [ab]n, 
one easily concludes that V* has an additive monodromy I(V) when passing around 
Cq counterclockwise. Moreover, as V G [0; 1] everywhere in A^, the boundary values 
of V* increases when going counterclockwise along Cn, as well as along dfl (recall that 
V* satisfies Neumann boundary conditions on and [afo]n). 

Let an additive constant in definition of V* be chosen so that V* = on dfl\ [ab]n 
(on some sheet of Aq). Then, there exists a non-self-intersecting nearest-neighbor path 
7 running from Cn to dfl \ [ab]n in A^ which separates non- negative (to the left of 7) 
and non-positive (to the right of 7) values of V*. We cut A^ along 7 and choose a 
branch of V* in A Q so that 

V* ^ at faces touching 7^ ht , V* ^ I(V) at faces touching 7^, 

V* = at faces touching [da] n , V* = I(V) at faces touching [bd] n 

(recall that V* satisfies Neumann boundary conditions on [aft]n). Putting on dual edges 
(//') = (vv'y of Al a discrete metric g ff , := \V*(f) - V*(f)\ = w vv ,\V(v') - V(v)\, 
one obtains the following estimate for the dual discrete length L* (see Remark 6.5) 
between opposite sides 7 nght U [da]n and [bd]n U 7 left of A Q : 

l* > viv)? > \mi = 1I(V) 

^ E (vv r )inAl ^v>\V(v>)-V(vW - 2I(V) 2 1 ) 

(the constant 2 is a big overkill, since each edge of Aq except 7 is counted once in A n , 
and only those constituting 7 are counted twice). Therefore, 

L A ,(Cn; [o6]n) = (L*)- 1 ^ 2[I(V)]- 1 = 2L An (C n ; [ab] a ) . 



(ii) Let d G dfl \ [ab]n be a boundary vertex where V attains its maximum on dfl 
(recall that V = on [ab] n ). If V(d) < then the metric g vv > := \V(v')-V(v)\ 

(which is extremal for the family (A^; -H- [afe]n), see Proposition 6.4) provides an 
estimate 



L An (Cn;dfl) >j^ = ^ 2 L As2 (Cn; [ab] Q ). 

If V(d) ^ 1 — g" 1 , let 7 denote a nearest-neighbor path running from d to C^ such 
that V ^ 1 — q^ 1 along this path (7 exists due to the maximum principle). Then, the 
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same metric as above (we assign zero weights to all edges constituting y eft ? ^ ri g ht and 
corresponding boundary ones) gives 

La^tS? UC.U 7**; [ab] a ) > = (l-g" 1 ) 2 L An (C n ; [ab] n ). □ 

Combining estimates given above, we are now able to prove a uniform double- 
sided estimate relating the logarithm of discrete harmonic measure ooq(u; [ab]n) in a 
simply connected fl and discrete extremal length L Af! (Cf}; [ab]n) in the annulus Aq(u). 

Theorem 7.8. Let a simply connected discrete domain Q and u G Int Q be such that 
Aq(u) is doubly connected, and [ab]n C dQ. Then, 

log(l + (w n («; [ab]^)- 1 ) x L A(l(u) (Cn(u); [06] n ), (7.5) 
with constants independent offl, u and \ab\o_. 

Proof. Let L := L Ao (Cn; [ab]n) and u := ujq(u; [ab]si). Corollary 7.4, Theorem 7.1 and 
Proposition 7.6 provide us the following diagram (for some proper discrete cross-cuts 
7 which can be different for lower and upper bounds): 

const .Z AJ r(C n ;[a6] n ) < u < Z A , (7^ U C^ U 7 ^ ght ; [ab] n ) 

t t 
2L > L A ,(C n ;Mn) L A , ( 7 >f U C Q U 7 " d ght ; Mn) > | L 

(the last inequality holds true, if L ^ Ao, where Ao is some absolute constant: recall 
that L As2 (Co;9n) x 1 for all Q and u). Above, u \" means a double-sided estimate of 
Z a t via L A 7 given by Theorem 7.1 (note that it is inverse monotone: an upper bound 
for L A 7 gives a lower bound for Z A 7 and vice versa). 

In particular, if L ^ Ao, condition (7.1) holds for both (right, and therefore, left) 
columns. Thus, in this case, one can replace both by (7.2), arriving at logo; x — L. 
If L < Ao, than the left column gives oj ^ const and both sides of (7.5) are uniformly 
comparable to 1 (note that L is uniformly bounded below by L Af2 (Cn; dfl) x 1). □ 

Corollary 7.9. Let Q be a simply connected domain, u G Intf2 ; and [ab]n G dfl. 
Denote oo disc := oo n (u; [ab] n ) and co cont := Uq(u; [ab] 1 ^). Then, 

log(l+a)^)xlog(l+0 

with some absolute (i.e., independent of fl,u,a,b) constants. 

Proof. First, let us assume that Aq(u) is doubly connected, so Q and u fit the setup of 
Theorem 7.8. Let Ldi SC : = L An (C^; [ab]n) and L cont := L a c(Cq; [ab]^) be its continuous 
counterpart. Recall that Ldi SC x L cont due to Proposition 6.2 (and Remark 7.5). Then, 

Mi+^disc) X L discr X L cont X log(l+UJ-^ nt ), 

where the last estimate is an easy corollary of the classical estimates for harmonic 
measure via extremal lengths (e.g., see [GM05, Theorem 5.2]). 

If Aq(u) is not doubly connected, then at least one of vertices v G dB n (u) belongs 
to a face touching dfl. In particular, this gives r Vint x dn(vi n t) x dn(u). Thus, the 
number of vertices #B^(m) is uniformly bounded (see Remark 2.2(h)). Further, if 
at least one of vertices of <9B^(-u) shares a face with [a6]n, then w disc ^ const and 
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w CO nt ^ const as well (indeed, the Brownian motion started at u has a positive chance 
to hit [ab]c traveling nearby of the corresponding finite length lattice path). 

Thus, without loss of generality, we may assume that both Wdisc and u cont are 
uniformly bounded away of 1 and there exists a connected (and simply connected) 
component of Int Aq(u) whose boundary contains the whole arc [ab]n. Let A' n denote 
this component of and Cq C <9Aq be the corresponding part of Cq slightly enlarged 
so that it includes two nearby boundary points of dQ (see Fig. 3B). Further, let L^ isc : = 
LAj J (Cf 2 ; [ab]a) and L' cont : = L A /c(Cq C ; [o&]H) denote its continuous counterpart. It is 
easy to see that one still has 

Wdiscr ~ Z n (w; [ab]n) x Z A ^(Cn n dA' n ; [ab] n ) x Za^(C^; [ab] n ) 

(the proof of Lemma 7.3 works well without any changes, and replacing H dA' n 
by costs no more than an absolute multiplicative constant). Applying (7.2) and 
Proposition 6.2, one obtains 

log w disc x -Lj iscr x -U coat x log u cont , 

(to prove the last estimate, e.g., draw a circle c u C Q c of radius \r u x dn(u) around u, 
then - log u cont x L n c(c u ; [ab]^) x L^ ont ). □ 

A. Appendix 

In this Appendix we show how to derive several statements of discrete potential 
theory used in the main text from two Assumptions (S) and (T) made in Section 2.3. 

A.l. Assumption (S), elliptic Harnack principle, and weak Beurling esti- 
mate. We start with a simple Lemma A.l which is an immediate corollary of (S), and 
then use it to derive (uniform) elliptic Harnack and weak Beurling estimates. 

Lemma A.l. There exist two constants 5 ,c > such that the following holds true: 
if u G T, r ^ r u are such that r v ^ 5 r for all v G B^(-u) (cf. Remark 2.2), and 
w G T : |r ^ \w — u\ ^ |r 7 then the probability of the event that the random walk 
started at w makes a whole turn inside of the annulus B^(-u) \ BF, 2 (u) (and cross its 
own trajectory afterwards) is uniformly bounded below by c > 0. 

Proof. Indeed, there exists a small constant tq > (depending on constants in assump- 
tions (a)-(c) in Sect. 2.1 only) such that applying (S) to a uniformly bounded number 
discs of radius r r and appropriate boundary arcs of angle n— r] step by step, one can 
"drive" the random walk started at w so that it makes the full turn staying inside of 
B£(u) \ B T r/2 (u) (cf. [CS11, Proof of Proposition 2.11]). □ 

Lemma A. 2 (elliptic Harnack principle). Let H : Q — > [0, +oo) be a nonnegative 
discrete harmonic function defined in Q. Then, for any u G Intf2 7 

H{y) x H{y') uniformly for all v, v' G Int B^^^-u), 

where dn(u) denotes the (Euclidean) distance from u to dVL. 

Proof. Note that, without loss of generality, we can assume that r v ^ 5 dn(w) for all 
v G B^(-u) (see Remark 2.2, otherwise the number of vertices in B^,Ju) is uniformly 
bounded, hence all values of H in this disc are uniformly comparable to each other, 
as all vertices can be joined by a uniformly bounded number of edges). Further, by 
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the maximum principle, there exists a nearest-neighbor path r y v running from v to dfl 
such that H > H(v) along this path. Then, Lemma A.l implies H(w) ^ c H(u) for 
all w : |dn(w) ^ \w — u\ ^ |d^(u) (as the random walk traveling around the annulus 
should cross 7^ at some point). Then, the maximum principle gives H(v') ^ CqH(v) 
for all v' G Int B^, w 2 (tx). The inverse estimate is similar. □ 

Remark A. 3. Using the same argument, one can prove the following one-sided version 
of Harnack's estimate near the boundary of a discrete domain fl which was used in the 
proof of Lemma 3.2 (recall that Bp(x) denotes the r-neighborhood of x in fl): 

Let x G dfl and r ^ r x be such that r v ^ 5 r for all i> G B^(x), and 
u> G B^(x) be such that |r ^ |w — tt| ^ |r. If if is a nonnegative 
discrete harmonic function in Bj?(x) which vanishes on dfl fl 9B^(x), 
then if(u) ^ c^Hiw) for all u G IntB^ /2 (x). 

Indeed, due to the maximum principle, one has H ^ H(v) along some nearest-neighbor 
path 7„ running from v to the outer boundary of B^(x) (this path cannot end on 
dfl fl <9B^(:r) as if vanishes there). Thus, Lemma A.l implies H(w) ^ CqH{v) since 
the random walk started at w has a positive chance to hit 7^ before dfl. 

Lemma A. 4 (weak Beurling estimate). There exists an absolute constant (3 > 
such that, for any simply connected discrete domain fl, u G Int fl and E C dfl, the 
following is fulfilled: 



ojq(u; E) ^ const 
where 



dn(u) 



distn(u;E) 



and un(u; E) ^ const 



diamE 
distn(w; E) 



distn(w; E) :— inf{r : u and E are connected in fin B^(u)}. 
Above we set diami? := r x , if E = {x} consists of a single vertex. 

Proof. The proof mimics the proof of Proposition 2.11 in [CS11]. Using Lemma A.l, it 
is easy to conclude that the random walk started at u has a chance Co > to hit dfl \ E 
before crossing the annulus B£ fcdn(u) \ B^ fc _ ldfi(u) for any k < log 2 [distn(w; E)/d n (u)] (if 

some annulus contains a vertex which local scale size is comparable to 2 k df l (u), then 
there are only a uniformly bounded number of vertices in B^^^, and so the random 
walk can hit dfl \ E by a uniformly bounded number of steps). This gives the first 
estimate with (3 := — log 2 (l — cq), and the second follows from the same arguments 
(applied to the random walk started on E and running to u). □ 

A. 2. Assumption (T) and the pointwise estimate for Green's function. For a 

simply connected domain fl, u G Int F and r > 0, denote 

B ( r\u):=B r rdn(u) (u) and C%>(u) := {x int G T : (x; (x int x)) G dB%>}. 
Also, let 

G^\v;u) := G (r)(v;u) 
denotes Green's functions in the discrete disc of radius rdn(u) around u. 
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Lemma A. 5. Let Q be a simply connected discrete domain and u G Intfi. Then, 

(i) for any v G Int B^\u), one has 

G^\v;u)^Gn(v;u)^c^GS\v;u); 

(ii) for any given < r < 1, there exist two constants c(r),C(r) > (independent of 
Q and u) such that the following holds true: 

c(r) < G n (v, u) ^ C{r) for all v G c£\u). 

Proof. Note that, without loss of generality, we can assume that r v ^ 5dn(w) for 
all v G B[^(w), where 5 is chosen small enough (otherwise, there is only a uniformly 
bounded number of points in B^(-u), so all values of Green's functions G^\ Gq, G^ are 

uniformly comparable with those at v — u, and Gq\u; u) x Gq(u; u) x Gq\u; u)xl 
since the random walk started at u has a positive chance to hit boundary by a uniformly 
bounded number of steps before coming back to u next time). 

(i) The first estimate G$ ^ Gn is trivial as Green's function is monotone with respect 
to Q, so we need to prove that Gq ^ c 1 G'[^' ) . Let v G Cq\u) for some r ^ 1. Using 
Lemma A.l, it is easy to conclude that 

G n (v;u) ^ G^\v;u) + (l-c ) max^^ G n (v;u). 

Indeed, if the random walk started at v reaches dB^\u) before hitting dQ, then it has 

a positive chance c > to hit dQ before coming back to Cq\u) (in order to have a 
possibility to visit u and to contribute to Gq(v;u) again). Thus, for any r ^ 1, 

max , e c s M(«) G n( v > u ) ^ c o 1 max , e c<[V) u ) ■ 

(ii) To prove the upper bound in (A.l), note that, for any v G Cq\u), one has 
Gci( • ; u) ^ Gq(v ; u) along some nearest-neighbor path running from v to u. As above, 
applying Lemma A.l and the maximum principle, we get 

Gn{w;u) ^ c Gn{v;u) uniformly for w G Int B^^ 6 (u), (A.l) 

(2) 

which implies Gq(v; u) ^ const due to G n (w; u) ^ cqG^{w] u) and the uniform upper 
bound in Assumption (T). 

To prove the lower bound in (A.l), note that m r := m^^H, , Gq\v ; u) is 
uniformly bounded below, if r is small enough. Indeed, the maximum principle gives 

Gq\v; u) ^ m r outside of Cq\u) , 

G$ (v; u) ^ m r + G^ (v; u) inside G$ (u) . 

Hence, Assumption (T) applied for both discs B^\u) and B^(-u) implies 

const •(dn('u)) 2 ^ const • m 2 (dn(v)) 2 + const •(rdn(f)) 2 . 

Therefore, there exist a constant r > (independent of Q and u) such that 
G^\v;u) ^ const > at least at one point v G Cq°\u). Applying Harnack's principle 
(Lemma A. 2) several times (depending on r and r but not on Q and u), one gets the 
uniform lower bound for G^( ■ ; u) ^ • ; u) on C^(w) for any given 0<r < 1. □ 
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Remark A. 6. In the proof of Proposition 3.1, we used the following uniform estimate: 
^ r 2 , G n (v; w) x (d n (u)) 2 for any w G Int B n (u). 

u6lntB n («) 

The lower bound is an easy corollary of the lower bound in (A.l). The upper bound 

(2) 

follows from (i) and the upper bound in (A.l) applied for G n (v;w). 
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